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I.  Introduction 


As  technology  advances,  thennal  management  challenges  continue  to  increase  in 
every  area,  from  multi-million  dollar  war  fighters  to  home  computer  systems.  New 
technology  is  both  more  powerful  and  more  compact.  The  by-product  of  higher  power  in 
a  smaller  package  is  that  an  immense  amount  of  heat  over  a  small  area,  also  known  as 
heat  flux,  must  be  dissipated  to  keep  systems  running  safely.  Mudawar  (2000)  identifies 
two  heat  flux  ranges  relative  to  the  amount  of  heat  dissipation.  The  high-flux  range 
includes  heat  fluxes  on  the  order  of  1 02- 1 03  W/cnr  and  the  ultra-high- flux  range  includes 
requirements  of  10  -10'  W/cm“.  Depending  on  the  fluid  and  cooling  method,  heat 
transfer  coefficients  for  removing  this  heat  may  range  from  0.001  to  100  W/(cm“-K). 
Single-phase,  forced  convection  improves  on  natural  convection  by  forcing  a  fluid  over 
the  heated  surface,  thus  producing  higher  heat  transfer  coefficients.  Boiling  (i.e.,  phase 
change  cooling)  takes  advantage  of  phase  change  of  the  liquid,  using  the  latent  heat  of 
vaporization  to  further  improve  on  the  attainable  heat  transfer  coefficients  and  increase 
the  critical  heat  fluxes  (CHF’s)  in  cooling  systems.  Compounding  the  challenge  is  the 
issue  of  varying  gravity,  which  is  found  in  many  aircraft  and  space-based  applications. 
Every  cooling  scheme  is  affected  in  some  way  by  gravity,  or  a  lack  thereof,  and  finding 
ways  to  deal  with  changes  in  gravity  is  vitally  important  to  building  effective  thennal 
management  systems  for  such  platforms. 

Predominant  cooling  fluids  include  air,  water,  and  fluorochemical  liquids.  Air  is 
the  cheapest  and  most  readily  available  cooling  fluid,  but  has  relatively  poor  thermal 
transport  properties,  so  its  use  is  limited  to  lower  heat  flux  requirements.  Water  remains 
the  most  effective  cooling  fluid  because  of  its  superior  thermal  transport  properties. 
However,  electrical  conductivity  precludes  its  use  in  direct  contact  with  electrical 
components  and  systems.  Consequently,  dielectric  fluorochemical  liquids,  such  as  those 


1 


produced  by  3M,  are  necessary  for  use  with  current-carrying  components,  despite  the  fact 
that  their  thermal  properties  are  inferior  to  those  of  water  (Mudawar,  2000). 

Phase  change  cooling  can  exist  in  several  fonns,  or  cooling  schemes.  Pool 
boiling  may  be  used  in  thennosyphons,  which  passively  cool  by  relying  on  buoyancy 
forces  to  circulate  coolant.  The  device  to  be  cooled  is  submerged  in  a  dielectric  coolant, 
which  boils  as  the  device  dissipates  heat,  releasing  vapor  bubbles  which  are  driven  by 
buoyancy  forces  to  the  upper  portion  of  the  container.  Here  the  vapor  condenses  and  falls 
back  into  the  coolant  pool.  Thermosyphon  performance  can  be  enhanced  by  several 
techniques,  including  surface  enhancement,  subcooling  the  fluid,  and  increasing  coolant 
pressure.  While  thennosyphons  are  only  one  possibility,  pool  boiling  in  general  relies  on 
buoyancy  to  carry  vapor  bubbles  from  the  heater  surface.  Thus,  orientation  has  a 
significant  effect  on  perfonnance  for  any  pool  boiling  system.  In  addition,  gravity  affects 
performance  as  well.  One  example  of  this  is  a  thermosyphon  relying  on  gravity  to  drain 
the  condensate  back  to  the  fluid  pool.  Much  research  has  been  conducted  in  pool  boiling 
(Kim  et  al.,  2002;  Kim,  2003;  Straub,  2001),  including  the  effects  of  microgravity.  It  has 
been  observed  that  CHF  is  generally  reduced  in  microgravity  (Kim,  2003),  due  in  part  to 
the  variations  in  bubble  dynamics  resulting  from  the  lack  of  buoyancy  forces  (Straub, 
2001). 

Heat  pipes,  like  thennosyphons,  are  passive  cooling  devices.  A  detailed 
description  of  heat  pipes  is  given  by  Faghri  (1995).  A  heat  pipe  consists  of  three  primary 
sections:  evaporator,  condenser,  and  transport  (adiabatic)  section.  The  coolant  begins  in 
the  evaporator,  where  it  is  exposed  to  the  heat  source  via  the  pipe  wall/wick  structure, 
and  evaporates,  building  pressure  within  the  vapor.  A  pressure  gradient  drives  the  vapor 
through  the  transport  section  to  the  condenser.  The  condensation  of  the  vapor  releases 
the  latent  heat  of  vaporization  to  the  ultimate  heat  sink.  Capillary  pressure  then  pumps 
the  condensate  through  the  wick  structure  back  to  the  evaporator,  where  the  cycle  begins 
again.  Heat  pipes  exist  in  many  forms  and  have  a  wide  range  of  applications.  Micro  heat 
pipes  may  be  as  small  as  25  mm  in  length,  whereas  macro  heat  pipes  may  be  as  long  as 
100  m.  Some  may  not  contain  an  adiabatic  section,  but  all  have  an  evaporator  and 
condenser.  Much  research  has  been  perfonned  on  heat  pipes  with  respect  to  the  effects  of 
acceleration.  Kaya  (2003)  used  a  centrifuge  table  to  produce  elevated  gravity  conditions, 
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and  detennined  that  the  accelerating  forces  at  high  spin  rates  can  cause  a  drying  out  of  the 
wick,  since  the  acceleration  overcomes  the  pumping  capacity  of  the  wick.  Gu  et  al. 
(2004)  studied  aluminum  pulsating  heat  pipes  (PHP)  in  microgravity.  They  determined 
that  PHP’s  are  capable  of  operating  successfully  in  microgravity,  and  may  be  suitable  for 
space-based  applications  using  lower-power  requirements. 

Active  cooling  systems  can  alleviate  some  of  the  issues  found  in  passive  systems 
and  improve  on  heat  transfer  characteristics.  They  are  more  suitable  for  high-  and  ultra- 
high-flux  applications.  Flow  boiling  sweeps  bubbles  off  of  heated  surfaces,  producing 
higher  heat  transfer  rates  and  higher  CHF  levels  (Kim,  2003).  As  in  pool  boiling, 
subcooling  the  liquid  helps  to  increase  CHF,  but  flow  boiling  has  the  added  ability  to 
vary  flow  velocity.  Increasing  flow  velocity  also  helps  to  raise  CHF,  in  addition  to 
reducing  the  significant  effects  of  flow  orientation  (Mudawar,  2000).  It  is  not  fully 
known  how  low  gravity  affects  flow  boiling,  as  power  requirements  and  large  transients 
have  prevented  much  research  in  this  area  thus  far. 

Jet  impingement  cooling  involves  directing  coolant  from  a  circular  or  rectangular 
orifice  onto  a  heated  surface.  Because  of  the  impact  momentum  involved,  it  is 
considered  an  aggressive  cooling  scheme,  and  is  not  feasible  in  cases  involving  delicate 
devices  or  components.  There  are  three  primary  formats  for  jet  impingement:  free  jet 
impingement,  where  the  liquid  is  introduced  into  a  gaseous  ambient,  submerged  jet 
impingement,  where  the  jet  is  issued  into  liquid  surroundings,  and  confined  jet 
impingement,  where  the  jet  is  confined  exclusively  between  the  orifice  and  heated 
surface.  Subcooling,  jet  diameter,  and  jet  velocity  may  all  increase  CHF  in  jet 
impingement,  but  jet  velocity  has  a  much  larger  effect  than  jet  width  (Mudawar,  2000). 
In  studying  jet  impingement  in  zero  gravity,  Labus  (1977)  defined  three  flow  regimes 
within  inviscid  flow  for  this  cooling  method:  surface  tension  flow,  transition  flow,  and 
inertial  flow,  where  each  regime  was  correlated  to  Weber  number  and  the  ratio  of  jet 
radius  to  disk  radius. 

Spray  cooling  has  been  used  for  many  applications  similar  to  those  using  jet 
impingement  cooling,  such  as  metal  quenching.  One  major  disadvantage  of  jet 
impingement  cooling  is  the  creation  of  large  temperature  gradients  within  the  device  due 
to  the  concentrated  heat  removal  in  the  vicinity  of  the  jet.  Spray  cooling  breaks  the 
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coolant  down  into  droplets,  which  individually  impinge  on  the  surface  and  provide  a 
more  unifonn  heat  removal.  Sprays  may  be  atomized,  where  air  (or  some  inert  gas)  is 
introduced  to  assist  in  the  breakdown  into  droplets.  Atomized  sprays  are  not  always 
practical,  however,  due  to  the  difficulty  of  removing  air  from  the  dielectric  cooling  liquid. 
Pressure  sprays  simply  supply  high  pressure  liquid  through  an  orifice,  and  do  not  use  air 
to  assist  in  the  droplet  formation  (Mudawar,  2000). 

Comparatively  little  research  has  been  done  in  the  area  of  spray  cooling  in 
microgravity.  Choi  and  Yao  (1987)  studied  heat  transfer  mechanisms  of  horizontally 
impacting  sprays.  They  observed  that,  as  the  liquid  mass  flux  increased,  so  did  the  heat 
transfer  capability.  Kato  et  al.  (1994)  used  water  and  CFC-113  as  test  fluids  in  spray 
cooling  systems  and  found  that  microgravity  produced  a  reduction  in  CHF  for  CFC-1 13, 
but  an  increase  in  CHF  for  water.  Yoshida  et  al.  (2001)  studied  both  copper  block  and 
thin- film  1TO  heaters  on  a  microgravity  platform,  and  found  that  gravity  dependency 
varied  with  spray  volume  flux  and  the  Weber  number  of  fluid  droplets  for  both  water  and 
FC-72. 

This  thesis  is  part  of  an  ongoing  joint  effort  between  the  Air  Force  Research 
Laboratory  (AFRL)  and  the  National  Aeronautics  and  Space  Administration  (NASA)  to 
study  spray  cooling  for  applications  in  both  reduced-  and  elevated-gravity  conditions. 
Previous  research  (Kato,  1994;  Yoshida,  2001)  indicates  spray  cooling  is  useful  to 
effectively  deal  with  waste  heat  generated  in  space-based  platforms,  such  as  the 
International  Space  Station  or  shuttle  orbiter,  as  well  as  the  military  war  fighter,  where  in¬ 
flight  maneuvers  can  produce  both  high-  and  reduced-gravity  conditions  for  transient 
periods. 

The  long-range  goals  of  the  experiment  include  macrothermophysics, 
microthermophysics,  and  fluid  management.  Macro-scale  studies  include  the  effect  of 
microgravity  on  the  critical  heat  flux  (CHF)  and  on  heat  transfer,  in  particular  the 
convective  heat  transfer  coefficient  due  to  the  spray.  Micro-scale  studies  include  film 
dynamics,  such  as  the  effects  of  both  reduced-  and  elevated-gravity  conditions  on  the 
liquid  film  formation,  surface  nucleation  and  vapor  bubble  formation,  and  droplet 
rebound  effects  (or  lack  thereof).  Where  Kato  (1994)  and  Yoshida  (2001)  studied  the 
effects  of  microgravity  on  a  single-nozzle  system,  and  Lin  and  Ponnappan  (2003)  built  a 
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closed-loop,  continuous-flow  system  for  use  in  terrestrial  gravity,  the  present  experiment 
aims  to  combine  these  efforts,  studying  a  continuous-flow,  closed-loop  spray  cooling 
system  in  micro-  and  elevated-gravity  conditions.  In  addition,  it  is  desired  in  the  long¬ 
term  to  integrate  a  multi-nozzle  array  into  the  system. 

The  objective  of  this  thesis  is  to  present  experimental  and  numerical  data  acquired 
for  a  spray  cooling  system  using  subcooled  liquid  in  a  terrestrial  gravity  environment, 
building  on  the  initial  design  and  testing  of  the  system,  as  discussed  in  Baysinger  et  al. 
(2004).  Baseline  tests  were  done  in  preparation  for  experiments  to  be  conducted  aboard 
the  NASA  KC-135  reduced  gravity  aircraft.  A  steady  state  numerical  model  was  built 
and  tested  using  ANSYS  finite  element  software,  and  a  parametric  study  was  performed 
to  determine  effects  of  volumetric  flow  rate,  heat  transfer  rate,  and  orientation  with 
respect  to  gravity  on  the  heat  transfer  within  the  spray  cooling  system.  The  numerical 
model  was  compared  to  the  experimental  data,  and  was  used  to  determine  the  heat 
transfer  coefficient  due  to  the  spray  and  to  predict  the  heat  transfer  due  to  conduction 
within  the  glass  heater  pedestal  assembly. 
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II.  Experimental  Design  and  Testing 


The  objective  of  the  experiment  was  to  determine  the  heat  transfer  coefficient  due 
to  spray  cooling  on  the  top  surface  of  a  glass  heater  pedestal  assembly.  An  ITO  heater 
layer  was  spray  deposited  onto  the  top  surface  of  the  glass,  which  was  used  to  provide 
ohmic  heating.  Thermocouples  were  embedded  into  several  layers  of  the  glass  to  provide 
temperatures  within  the  pedestal.  These  temperatures  were  then  used  in  conjunction  with 
the  numerical  model,  described  in  the  Numerical  Modeling  section  of  this  thesis,  to 
determine  the  heat  transfer  coefficient  of  the  spray.  While  this  thesis  focuses  on  the 
baseline  results  of  the  experiment  in  Earth’s  gravitational  field,  eventually  the  experiment 
will  be  repeated  aboard  the  NASA  KC-135  reduced  gravity  aircraft.  Microgravity,  lunar 
gravity,  Martian  gravity  and  high-gravity  (approximately  1.8  g)  conditions  will  be 
produced  by  following  a  parabolic  flight  trajectory.  Each  flight  will  consist  of  about  40 
parabolas,  split  into  sequences  of  10  parabolas,  with  a  turn  after  each  set  of  10.  The 
experimental  series  will  be  comprised  of  four  flights,  for  a  total  of  approximately  160 
parabolas.  Figure  1  shows  a  typical  flight  profile  with  turns  before  and  after  a  series  of 
parabolas.  This  figure  also  shows  the  variation  in  ambient  cabin  temperature  which  is 
consistent  with  the  variation  in  aircraft  altitude  during  maneuvering. 

The  experimental  apparatus  consisted  of  two  primary  components  (Baysinger  et 
ah,  2003).  The  first  component  was  a  spray  test  chamber  containing  two  glass  heater 
pedestals,  two  spray  nozzles,  a  primary  condenser  surface,  and  two  sump  configurations 
to  collect  the  liquid  and  condensate.  The  second  component  was  the  flow  loop  system 
consisting  of  two  flow  loops:  1)  an  FC-72  loop  to  manage  the  working  fluid  and  2)  a 
water  loop  to  cool  the  spray  test  chamber.  A  photograph  of  the  overall  experimental 
apparatus  is  shown  in  Figure  2. 

Figure  3  shows  the  spray  test  chamber,  which  was  a  cylinder  having  interior 
dimensions  20.2  cm  (7.97  in)  in  diameter  and  15.2  cm  (6  in)  in  length.  It  was  made  of 
standard  vacuum  components  to  provide  a  leak-tight  environment  for  the  spray-cooling 
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experiment.  Key  components  included  opposing  nozzles  to  allow  for  simultaneous  testing 
of  different  heater  orientations  with  respect  to  gravity,  a  cooling  shell  surrounding  the  test 
chamber,  and  a  wick  structure  surrounding  the  test  area  to  collect  and  direct  excess  fluid 
and  condensate. 

Each  glass  heater  pedestal  assembly  was  cylindrical  with  a  diameter  of  0.016  m, 
and  consisted  of  five  layers,  as  shown  in  Figure  4  and  Figure  5.  The  first  layer  was  a  thin- 
film,  transparent,  indium-tin  oxide  (ITO)  heater,  which  was  sputter-deposited  onto  the 
0.5-mm-thick  fused  silica  substrate  (Corning  AMFCD  1737  glass).  This  combination  was 
chosen  to  match  the  coefficient  of  thermal  expansion  as  closely  as  possible  between  the 
ITO  layer  and  the  glass  layer  to  prevent  delamination.  The  heater/substrate  was  mounted 
onto  a  1-mm-thick  wafer  of  ZN-7  Borosilicate  glass  using  optical  cement  (Summers  SK- 
9).  This  subassembly  was  then  adhered  to  an  identical  1-mm-thick  wafer  of  Borosilicate 
glass,  which  was  finally  cemented  to  the  main  glass  pedestal,  which  was  also  made  of  the 
ZN-7  glass.  The  purpose  of  the  two  intermediate  layers  was  to  allow  for  grooves  to  place 
thennocouples,  as  shown  in  Figure  6.  The  heater  pedestal  assembly,  as  fabricated,  also 
allowed  for  an  optical  path  for  visual  observation  of  the  heater  surface.  A  preliminary 
analysis  was  used  to  predict  temperatures  throughout  the  pedestal  in  order  to  determine 
the  best  placement  for  the  thennocouples. 

Various  nozzles  were  evaluated  using  a  two-axis  phase  Doppler  anemometer 
(Dantec)  to  detennine  spray  unifonnity  and  droplet  size  distribution.  Typical  droplet 
velocity,  volumetric  flux,  and  mean  droplet  diameter  distributions  for  the  selected  nozzles 
are  shown  in  Figure  7,  Figure  8,  and  Figure  9,  respectively.  The  distance  from  the  nozzle 
to  the  heater  surface  was  set  to  approximately  1  cm.  This  distance  was  based  on  a  spray 
cone  angle  of  60°  such  that  the  spray  contacted  only  the  heater  surface  with  minimal 
overspray.  The  motion  of  any  excess  fluid,  which  would  ordinarily  drain  to  the  bottom  of 
the  chamber/pedestal  in  terrestrial  conditions,  is  unpredictable  in  microgravity,  possibly 
floating  in  the  chamber  or  adhering  to  structures  due  to  surface  tension.  Regardless,  the 
test  fluid  needed  to  be  effectively  removed  from  the  chamber  for  recirculation.  This  was 
accomplished  by  placing  a  sump  configuration  at  each  nozzle.  The  sumps  served  to 
remove  excess  fluid  from  the  chamber  and  spray  surface  through  an  annulus  between  the 
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heater  pedestal  and  sump  wall.  Figure  10  shows  the  sump  and  the  heater  pedestal 
together. 

Vapor  generated  at  the  heater  surface  was  exhausted  out  past  the  nozzle  into  the 
test  chamber  interior  or  entrained  with  the  excess  liquid  and  carried  through  the  sump 
annulus.  Excess  vapor  and  liquid  in  the  chamber  were  controlled  by  an  annular  wick 
structure,  both  to  aid  in  condensation  and  to  contain  the  condensate.  The  wick  structure 
was  formed  by  lining  the  internal  chamber  wall  with  one  layer  each  of  100  and  150 
stainless  steel  mesh,  with  the  150  mesh  being  nearest  the  wall.  An  additional  cylindrical 
screen  was  placed  in  the  chamber  interior  to  provide  for  an  annular  space  between  the 
screen  and  the  chamber  wall.  This  annular  space  between  the  chamber  wall  and  inner 
screen  was  filled  with  stainless  steel  wool  to  help  confine  the  unconstrained  liquid  during 
microgravity. 

The  test  chamber  was  affixed  to  a  hinge  enabling  it  to  be  rotated  so  that  the 
nozzles  could  be  oriented  in  either  a  vertical  or  horizontal  position.  This  allowed  greater 
flexibility  in  planning  and  executing  experiments.  For  the  experiments  described  in  this 
thesis,  the  chamber  was  positioned  such  that  the  nozzles  remained  in  the  vertical 
orientation.  A  pressure  transducer  was  installed  to  measure  the  vapor  pressure  inside  the 
chamber.  A  tri-axis  accelerometer  (Columbia  SA307HPTX)  with  an  uncertainty  of 
±0.03-g  was  mounted  onto  the  test  chamber  to  record  the  transient  variation  in 
acceleration  throughout  the  flight  profile. 

Figure  1 1  shows  a  schematic  diagram  of  the  fluid  loops.  Each  loop  contains  a 
combination  of  pumps,  flow  meters,  pressure  transducers,  valves,  pressure  switches, 
reservoirs,  liquid-air  heat  exchangers,  and  interconnecting  plumbing  to  monitor  and  direct 
the  fluid  flow.  Temperatures  throughout  the  fluid  loops  were  monitored  using  Type  E 
thennocouples  which  were  0.159  cm  (0.0625  in)  in  diameter.  Key  thennocouple  locations 
include  the  entrances  and  exits  of  the  test  chamber,  reservoirs,  and  heat  exchanger 
assemblies. 

A  positive  displacement  scavenging  pump  (Tuthill)  was  used  to  pump  subcooled 
FC-72  fluid  from  the  sumps  through  a  series  of  three  air/liquid  heat  exchangers  to  a 
reservoir.  Pressure  transducers  were  placed  on  either  side  of  the  pump,  and  the  fluid 
passed  through  a  40  pm  filter  (Swagelock)  before  reaching  the  reservoir.  The  reservoir 
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contained  either  stainless  steel  wool  or  a  vane  structure  to  reduce  fluid  sloshing  during 
the  flight  profde  and  served  as  both  a  volume  for  the  scavenging  pump  to  dump  excess 
working  fluid  and  a  pickup  reservoir  for  the  nozzle  pumps.  Following  the  reservoir,  the 
fluid  branched  off  into  two  identical  flow  paths,  one  for  each  spray  nozzle.  Each  branch 
consisted  of  a  positive  displacement  pump  followed  by  a  15  pm  fdter  (Swagelock),  a 
flow  meter  (Omega),  a  pressure  transducer,  and  a  pressure  switch  (United  Electric). 
Electrically-actuated  ball  valves  (Hoke)  were  placed  at  the  entrance  and  exit  of  each 
branch,  to  permit  isolation  of  the  test  chamber.  All  testing  was  done  with  subcooled  fluid 
at  ambient  temperature;  however,  a  length  of  copper  tubing  was  placed  in  each  branch  for 
future  installation  of  a  reheater  that  would  allow  for  testing  of  the  spray  at  varying 
degrees  of  subcooling  and  saturated  conditions. 

The  single-phase  water  loop  maintained  the  temperature  of  the  exterior  wall  of  the 
test  chamber  at  a  constant  value.  The  design  and  components  were  similar  to  those  in  the 
FC-72  loop.  A  positive-displacement  pump  moved  the  water  from  an  assembly  of  two 
liquid/air  heat  exchangers  to  a  reservoir  via  a  140  pm  filter,  pressure  relief  switch,  flow 
meter,  and  pressure  transducer.  As  with  the  FC-72  loop,  the  reservoir  contained  stainless 
steel  wool  to  reduce  fluid  sloshing. 

The  experiment  framework  was  constructed  from  “2020”  and  “1010”  sizes  of 
aluminum  extrusions  and  hardware  (80/20,  Inc.).  Shelves  were  constructed  of  0.635  cm 
(0.25  in)  thick  aluminum  and  installed  for  the  purpose  of  mounting  equipment.  The  entire 
frame  rested  on  a  1.27  cm  (0.5  in)  thick  aluminum  plate,  which  served  both  as  a 
mounting  plate  and  a  spill  containment  pan. 

Detailed  structural  analyses  and  testing  were  conducted  to  satisfy  NASA  safety 
requirements.  Structural  integrity  under  various  acceleration  field  strengths  was  evaluated 
for  flight  on  the  KC-135  aircraft.  NASA  guidelines  required  that  all  framework  and 
equipment  be  able  to  withstand  a  9-g  load  in  the  forward  direction,  3-g  in  the  aft,  6-g 
down,  2-g  lateral,  and  2-g  up.  In  addition,  a  factor  of  safety  of  two  was  required,  so  all 
structural  analysis  and  testing  were  perfonned  at  twice  the  given  requirements.  Much  of 
the  smaller  equipment  installed  on  the  experiment  was  pull-tested  using  a  load  cell,  but 
larger  items  and  the  frame  itself  were  analyzed  using  static  equilibrium  equations. 
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The  experiment  was  operated  and  monitored  via  a  control  panel  and  data 
acquisition  system.  The  data  acquisition  system  consisted  of  a  laptop  computer  mounted 
to  the  top  of  the  experimental  package,  coupled  with  a  data  acquisition/switch  unit  (HP 
34970A).  A  diagram  of  the  instrumentation  system  is  shown  in  Figure  12,  and  a 
photograph  of  the  control  panel  is  shown  in  Figure  13. 

Temperature  measurements  were  obtained  throughout  the  chamber,  glass  heater 
pedestal,  and  within  the  sump  annulus  using  sheathed  thennocouples.  Type  E 
thennocouples  (0.16  cm  diameter,  time  constant  <  4.0  s)  were  placed  in  the  chamber  to 
measure  chamber  wall  and  vapor  temperatures,  and  Type  E  thermocouples  (0.0254  cm 
diameter,  time  constant  <  0.5  s),  were  mounted  within  the  glass  heater  pedestal  assembly. 
For  additional  stability  and  accuracy  in  temperature  measurement,  the  thermocouples 
within  the  chamber  and  glass  heater  pedestal  assembly  were  referenced  to  a  zero  point 
reference  junction  (Hart  9101)  with  a  stability  of  ±0.005°C. 

Temperatures,  as  measured  by  thermocouples,  within  the  glass  heater  pedestal 
assembly  were  used  to  determine  the  heat  transfer  coefficient  due  to  the  liquid  spraying 
onto  the  top  of  the  pedestals.  In  addition,  the  ambient  temperature  was  monitored  using  a 
Type  E  thennocouple  and  a  thermistor  (YSI  44006)  mounted  on  the  frame  of  the 
experiment.  Thermocouples  mounted  within  the  chamber  and  glass  heater  pedestals  were 
calibrated  using  an  RTD  temperature  probe  (Hart  Scientific  Model  5628)  and  temperature 
bath  (Hart  Scientific  Model  7320)  in  the  temperature  range  from  35°  to  120°C.  Resulting 
thennocouple  calibrations  yielded  an  uncertainty  of  ±0.02°C.  A  typical  thermocouple 
calibration  curve  is  shown  in  Figure  14.  Two  additional  Type  E  thermocouples  were 
placed  at  the  edge  of  the  heater,  to  measure  the  liquid  temperature  exiting  off  of  the 
heater  surface,  and  in  the  sump  annulus,  approximately  1 .27  cm  below  the  heater  surface. 

Heater  power  was  monitored  using  voltage  measurements  with  an  uncertainty  of 
±0.05%  across  the  heater  and  a  precision  resistor  with  a  resistance  of  0.10  Q  ±  0.02%, 
resulting  in  a  heater  power  measurement  error  ±0.07%.  Visual  documentation  was 
acquired  through  the  use  of  two  video  cameras  (Sony  DCR-TRV33)  mounted  to  the 
experimental  frame.  One  camera  recorded  a  full  view  of  the  test  chamber,  observing  the 
behavior  of  the  fluid  therein.  The  second  camera  was  positioned  to  record  a  field  of  view 
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of  only  the  nozzle  being  tested.  When  testing  switched  from  one  nozzle  to  the  other,  the 
camera  position  was  adjusted  to  record  the  appropriate  nozzle. 

At  the  beginning  of  each  experiment,  the  data  acquisition  program  was  initiated 
via  a  laptop  computer,  and  all  cameras,  timers,  etc.,  were  synchronized  and  started.  The 
data  acquisition  system  monitored  and  recorded  the  flow  rate,  acceleration,  input  heat, 
and  temperatures  for  the  duration  of  the  experiment.  Temperatures  were  measured  within 
the  glass  heater  pedestal,  at  the  inlet  of  the  nozzle,  and  in  the  liquid  film  exiting  the  heater 
surface.  During  the  experiments,  the  volumetric  flow  rate  and  heater  power  were  varied 
for  values  of  V  =  4,  5,  6,  8,  and  10  gal/hr  and  Q  =  20,  30,  and  50  W,  respectively. 
Experiments  were  conducted  for  a  fixed  flow  rate  while  varying  the  heater  power.  When 
the  heater  interface  temperature  reached  7j  =  100°C,  the  flow  rate  was  adjusted  and  the 
heater  power  varied  for  the  new  flow  rate  setting. 
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T  (°C) 


t(s) 

Figure  1:  Typical  flight  profile  with  ambient  temperature  and  acceleration  versus  time. 
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X- Accel  (g) 


Figure  3:  Solidworks  model  of  the  spray  test  chamber. 
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(b) 

Figure  4:  Details  of  the  glass  heater  pedestal  assembly:  (a)  Glass  parts  used  in  assembling 
the  heater  pedestal;  (b)  Assembled  heater  pedestal  with  embedded  thermocouples  and 
copper  wires  for  ohmic  heating. 
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(b) 


Figure  5:  Details  of  the  glass  heater  pedestal  assembly,  continued:  (a)  Side  view  of  the 
top  of  the  assembled  heater  pedestal;  (b)  Top  view  of  the  assembled  heater  pedestal  prior 
to  application  of  the  electrically  conductive  silver  epoxy. 
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Section  at  top  surface 


Section  at  1  mm  depth 


Section  at  2  mm  depth 


Section  at  50  mm  depth 


Figure  6:  Groove  locations  in  the  glass  wafers  for  thermocouple  placement  (dimensions 
in  mm). 
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[mm] 


Normal  Component  of  Velocity  [m/s] 
Spraying  Systems  FullJet  1  Nozzle 
9.0  gph,  40  psig,  15  mm 


X  [mm] 

Figure  7:  Normal  component  of  velocity  for  the  nozzle  (Spraying  Systems  FullJet  1, 
V  =  9.0  gal/hr,  p  =  40  psig,  h\  =  0.015  m  from  orifice). 
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Flux  vs  Position 

Across  the  Center  of  the  Spray 
for  a  Spraying  Systems  Fulljet  1  Nozzle 


Figure  8:  Volumetric  flux  versus  position  across  center  of  spray  (Spraying  Systems 
FullJetl  nozzle,  h}  =  0.017  m  from  the  orifice). 
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Droplet  Diameter  vs  Position 
Across  the  Center  of  the  Spray 
for  a  Spraying  Systems  Fulljet  1  Nozzle 


Figure  9:  Mean  droplet  diameter  versus  position  across  center  of  spray  (Spraying  Systems 
FullJetl  nozzle,  h}  =  0.017  m  from  orifice). 
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Figure  10:  Solidworks  model  and  photographs  of  the  heater  pedestal  assembly  and  sump 
configuration. 
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Nozzle  Flow  Meter 
&  Inlet  Fluid  Temperature 


Figure  12:  Schematic  diagram  of  the  instrumentation. 
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Figure  13:  Photograph  of  the  control  panel. 
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RTD  Temp  (C) 


III.  Mathematical  Model  of  the  Glass  Heater  Pedestal 
Assembly 


The  objective  of  this  section  is  to  mathematically  describe  the  glass  heater 
pedestal  assembly,  which  was  fabricated  as  a  series  of  five  layers  as  shown  in  Figure  15. 
The  ITO  layer  used  for  ohmic  heating  was  sputter-deposited  onto  the  first  layer  of  glass 
substrate.  Two  additional  thicker  layers  of  glass  substrate  were  machined  with  grooves  to 
accept  thermocouples  to  measure  temperatures  within  the  assembled  glass  heater  pedestal 
assembly.  The  three  layers  of  glass  were  mounted  onto  the  glass  pedestal  base  using 
optical  cement.  The  entire  assembly  was  placed  into  the  spray  chamber  under  the  spray 
nozzle.  The  liquid  spray  impinged  onto  the  top  surface  of  the  ITO  heater.  The  liquid 
droplets  were  able  to  either  evaporate  or  flow  off  of  the  surface  due  to  the  combined 
effects  of  gravity  and  the  momentum  of  the  spray.  While  part  of  the  spray  was  able  to 
rebound  off  of  the  heater  surface,  most  of  the  liquid  ran  off  the  top  surface  and  down  the 
pedestal,  either  in  rivulets  or  as  a  liquid  film.  For  the  ITO  heater  layer,  assuming  that  the 
thennal  properties  are  constant,  the  transient  conservation  of  energy  equation  is 
(Rohsenow  et  al.,  1985) 


PlC 


1W1 


8T{ 

dt 


=  k, 


"l  0 

(rd_TC 

r  dr 

V  dr  , 

1  8%  82TX 
+  ^r — ^  +  - 
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+  q  i 


where  q{  is  the  internal  heat  generation  per  unit  volume  due  to  ohmic  heating.  On  the  top 
surface  of  the  ITO  heater,  the  heat  transfer  coefficient  is  assumed  to  vary  along  the  radial 
and  azimuthal  directions  due  to  the  fact  that  the  spray  may  not  be  unifonn.  Also,  the  heat 
transfer  coefficient  is  potentially  transient  if  the  liquid  flow  rate  of  the  spray  varies  with 
time: 


0  <r<Rp,  0  <  0  <  27i,  z  =  Hp,  t>  0: 


8T 

-kx-^  =  hx(Tm-Ttl) 
dz 
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where  hx  =  hx  (r,0,z  =  H  ,t)  is  the  heat  transfer  coefficient  at  the  top  surface,  Tw  is  the 

local  wall  temperature,  and  Tf\  is  the  local  film  temperature  of  the  liquid  flowing  over 
the  heater  surface.  On  the  bottom  face  of  the  ITO  heater,  the  temperature  and  heat  flux 
are  continuous  due  to  the  perfect  thermal  contact  of  the  heater  with  the  first  layer  of  glass 
substrate,  which  is  a  result  of  the  sputtering  process,  giving  a  boundary  condition  of 


dT  dT 

0  <  r  <  R  ,  0  <  0  <  2n,  z  =  H1+H3  +  H  4  +H5,  t  >  0  :  kx  — -  =  k2  — — 

dz  dz 


T  =  T 
‘  1  1  2 


Within  the  first  layer  of  glass  substrate,  the  energy  equation  is  given  as  follows,  again 
assuming  that  the  properties  of  the  glass  do  not  vary  with  temperature. 

52t n  a2 r 


dT2 

P2Cp,2  a  —  *2 

Ot 


"1  d 

(  dTA 

r  dr 

l  dr  J 

1  8%  d% 

+  — r - f-  +  - 


dz" 


The  first  layer  of  glass  substrate  is  adhered  to  the  second  thicker  layer  of  glass  with 
optical  cement.  Therefore,  the  boundary  condition  at  the  bottom  face  of  the  first  glass 
layer  can  be  modeled  using  the  contact  conductance  of  the  interface  (Ozisik,  1993): 

dT  dT 

0  <r<R,  0  <  9  <  271,  z  =  H,+H4+H5,  t>  0:  -  k,  =  hc2(T2 -T3)  =  -k3 — 

dz  '  dz 

where  /zc, 2  is  the  contact  conductance  at  this  interface.  For  the  second  layer  of  glass,  the 
conservation  of  energy  equation  is 


p3c 


dT 


P.3 


dt 


=  k. 


"1  d 

(  dTA 

r  dr 

y 

l  dr  ) 

1  d%  d% 

+  — r - f-  +  - 


r  dQz 


dz < 


Again,  the  boundary  condition  on  the  lower  face  of  the  second  glass  layer  can  be  modeled 
using  the  contact  conductance. 


ZT  QT 

0  <r  <  R  ,  O<0<  271,  z  =  H4+H5,  t  >  0  :  -k3  — -  =  hc  3  (T3  ~T4)  =  -k4  — — 

dz  ’  dz 

For  the  third  layer  of  glass,  the  conservation  of  energy  equation  is 
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The  boundary  condition  on  the  lower  face  of  the  third  glass  layer  can  again  be  modeled 
using  the  contact  conductance. 
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dT  dT 

0  <  r  <  R  ,  0  <  0  <  27i,  z  =  H5,  t>  0:  -  k4  — 1  =  hc4  {T4  -  T5 )  =  -k5  — 5- 

dz  dz 

The  conservation  of  energy  equation  for  the  glass  base  is 
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The  lower  face  of  the  glass  base  is  assumed  to  be  perfectly  insulated,  due  to  its  length  and 
low  thermal  conductivity. 


dT 

0  <r<Rn,  O<0<  2n,  z  =  0,  t  >  0  :  — -  =  0 
p  dz 

Along  the  outer  radius  of  the  ITO  heater,  the  glass  substrate  layers  and  the  glass  base,  the 
heat  transfer  coefficient  varies  due  to  the  complicated  transient  behavior  of  the  liquid  as  it 
falls,  either  in  rivulets  or  as  a  thin  liquid  film.  Therefore,  the  boundary  conditions  along 
the  outer  radius  of  these  regions  are  as  follows: 

ITO  heater: 

r  =  Rp,  0  <  0  <  271,  H2+H,+H4+Hs<z<Hp, 

First  glass  layer: 

r  =  Rp,  0  <  0  <  271,  H3+H4+Hs  <z<Hp-Hl, 

Second  glass  layer: 

r  =  Rp,  0  <  0  <  271,  H4+Hs  <z<H,+H4+H5, 

Third  glass  layer: 

r  =  Rp ,  0  <  0  <  27i,  Hs<z<H4  +H5,  t>  0:  -  k4^  =  h2(Tw -Tn) 

dr 

Glass  base: 

dT 

r  =  R  0  <  0  <  271,  0  <z<H5,  t>  0:  -kx—L  =  h2(Tyi-Ti2) 

dr 


dT 

t>0:  -kl—L  =  h2(Tw-T{2) 
dr 


dT 

t>  0:  ~k2^  =  h2(Tw-T{2) 
dr 


£\rn 

t>  0:  -k3—^  =  h2(Tw-T{2) 
dr 


dT. 
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where  h2  =  h2  (r  =  R  ,Q,z,t )  is  the  heat  transfer  coefficient  along  the  outer  radius  of  the 
glass  heater  pedestal  assembly,  Tw=Tw(Q,z,t)  is  the  local  wall  temperature  and 
7f2  =  T(1  (0,  z,t)  is  the  local  film  temperature  of  the  liquid  falling  down  the  side. 

The  initial  conditions  for  the  above  conservation  of  energy  equations  are  that  the 
temperature  throughout  the  entire  domain  is  assumed  to  be  constant: 

0  <r<Rp,  0  <  0  <  271,  0  <  z  <  Hp,  t  =  0:  Tx  =  T2  =  T3  =  T4  =  T5  =  T0 

where  T0  is  the  initial  temperature. 

For  the  present  analysis,  the  following  assumptions  have  been  made: 

1 .  Steady  state  conditions  prevail. 

2.  The  thermophysical  properties  of  the  first  glass  layer  correspond  to  the  Corning 
glass,  whereas  the  second  and  third  glass  layers  and  the  glass  base  are  made  of 
borosilicate  glass.  The  properties  of  the  ITO  heater  are  also  specified. 

k\  =km,  Pj  =  P  iTO  ’  Cp,l  =  Cp,ITO 
^2  =  keg  ’  P  2  =  P  eg  5  Cp,2  =  Cp,cg 

k}  -  k4  -  k5  -  kbg,  P3  =P4  =P5  =Pbg5  Cp.3  “CP-4  “CP,5  “Cp,bg 

3.  The  heat  transfer  coefficient  on  the  top  surface  of  the  ITO  heater  is  constant. 
hx  (r,  0,  z  =  H  ,t)=  constant 

4.  The  heat  transfer  coefficient  on  the  outer  radius  of  the  glass  heater  pedestal 
assembly  is  constant. 

h2  (r  =  Rp ,  0,  z,  t )  =  constant 

5.  The  contact  conductances  at  the  interfaces  of  the  glass  layers  are  infinite,  which 
means  that  both  the  temperatures  and  heat  fluxes  at  these  interfaces  are 
continuous. 

6.  The  temperatures  used  to  define  the  heat  transfer  coefficients  are  obtained  from 
the  experimental  data.  The  freestream  temperature  for  the  spray  heat  transfer 
coefficient  is  the  average  of  the  fluid  temperature  entering  the  nozzle  and  the 
upper  film  temperature. 

^f,l  =  ^  (^noz  +  Tur ) 
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7.  The  freestream  temperature  for  the  heat  transfer  coefficient  on  the  side  of  the 
heater  pedestal  assembly  is  the  average  of  the  upper  and  lower  film  temperatures. 

Ta=\(T^+Ttt) 
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Figure  15:  Schematic  diagram  of  the  glass  heater  pedestal  assembly. 
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IV.  Numerical  Model  of  the  Glass  Heater  Pedestal  Assembly 


The  objective  of  this  section  is  twofold:  First,  a  preliminary  analysis  will  be 
described  in  which  transient  temperature  distributions  within  a  right-circular  cylinder 
were  determined  both  analytically  and  numerically.  This  study  was  performed  in  order  to 
become  confident  with  the  computer  program  used  in  this  analysis  (ANSYS  5.4  finite 
element  software).  Secondly,  time-step  and  grid  independence  studies  will  be  outlined  to 
show  that  the  results  of  the  numerical  heat  transfer  analysis  for  the  glass  heater  pedestal 
assembly  are  grid  independent. 

Preliminary  Analysis 

Prior  to  modeling  of  the  glass  heater  pedestal  assembly,  the  results  of  an 
analytical,  two-dimensional,  transient  heat  conduction  model  of  a  right  circular  cylinder 
were  compared  to  those  of  a  numerical  finite  element  model  using  the  software  package 
ANSYS  5.4.  In  this  way,  expertise  was  developed  in  the  modeling  of  heat  conduction 
problems,  which  could  then  be  used  to  model  the  glass  heater  pedestal  assembly. 
Beginning  with  the  energy  equation  in  cylindrical  coordinates,  appropriate  assumptions 
were  made  and  boundary  and  initial  conditions  were  applied  to  determine  the  solution 
analytically.  The  corresponding  numerical  model  was  then  developed,  and  transient 
temperature  profiles  within  the  right  circular  cylinder  were  obtained.  The  results  of  the 
analytical  model  and  those  of  the  numerical  model  were  compared,  and  it  was  found  that 
the  two  models  were  in  agreement. 

A  schematic  diagram  of  the  analytical  model  of  the  transient  heat  conduction 
within  a  right  circular  cylinder  of  finite  length  is  shown  in  Figure  16(a).  The  energy 
equation  for  this  case  is  given  by  Rohsenow  (1985): 

pcn  —  =  V-kVT  +  q"'  (4.1) 

v  dt 
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It  is  assumed  that  all  of  the  thermal  properties  are  constant,  there  is  no  internal  heat 
generation,  and  that  the  solution  domain  is  axisymmetric,  as  shown  in  Figure  16(h). 
Under  these  assumptions,  the  energy  equation  reduces  to 

d2T  1  dT  d2T  1  dT 

— r  + - +  ^  = -  (4-2) 

dr2  r  dr  dz 2  a  dt 

On  the  top  face  of  the  cylinder,  a  constant  heat  flux  is  applied. 

0  <r<R,  z  =  H,  t>  0:  —  =  (4.3) 

dz  k 

On  the  bottom  face,  the  temperature  is  zero. 


0  <r<R,  z  =  0,  t>  0:  T  =  OK  (4.4) 

At  the  centerline,  the  symmetry  condition  is  imposed. 

r  =  0,  0  <z<H,  t>  0:  —  =  0  (4.5) 

dr 

At  the  outer  radius  of  the  cylinder,  the  temperature  is  zero. 

r  =  R ,  0  <z<H,  t>  0:  T  =  OK  (4.6) 

At  the  initial  time,  the  temperature  throughout  the  cylinder  is  constant. 

0  <r<R,  0  <  z  <  H,  t  =  0:  T  =  T  (4.7) 

The  energy  equation  and  the  boundary  and  initial  conditions  are 
nondimensionalized  using  the  following  definitions: 


T 

T*  =  —, 
T 


*  z  *  r  ^  at 

z*  =  — ,  r*  =  — ,  t*  = T- 


The  conservation  of  energy  equation  becomes 

d2T  *  1  dT  *  d2T  *  dT  * 

- —  ~ ^ - 1 - —  = - 

Qr  *  r  *  Qr  *  Qz  *-  fa  * 

The  boundary  conditions  and  the  initial  condition  reduce  to 


0<r*<B,  z*  =  1,  t*  >  0 : 


(4.8a) 


0<r*<B,  z*  =  0,  >  0  :  T*  =  0 


(4.8b) 
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r*  —  0,  0  <  z*  <  1,  t*  >  0  : 


sr* 

dr  * 


=  0 


r*  =  B,  0  <  z*  <  1,  t*  >  0 :  T*  =  0 
0  <r*<B,  0  <  z*  <  1,  f*  =  0 :  T*  =  1 


where  the  aspect  ratio  is  B  =  R!  H  and  the  dimensionless  heat  flux  is  given  by 


<h  H 
kTj 


(4.8c) 

(4.8d) 

(4.8e) 


(4.9) 


Due  to  the  linear  nature  of  the  above  dimensionless  problem,  it  can  be  divided 
into  a  transient  portion  and  a  steady-state  portion,  where  each  of  the  sub-problems 
contains  one  nonhomo geneous  condition  (Figure  17).  The  dimensionless  temperature  is 
defined  as  T*  =  Tx  +T2,  where  the  *’s  have  been  dropped  for  clarity.  The  steady-state 
problem  is  given  by  the  following  equations: 


a%  + 1  sr,  +  02r,  _  Q 

dr  r  dr  dz2 

(4.10) 

dT 

0  <r<B,  z  =  1 :  -i  =  -g0 

dz 

(4.10a) 

0  <r<B,  z  =  0  :  Zj=0 

(4.10b) 

dT 

r  =  0,  0  <  z  <  1 :  ^  =  0 

dr 

(4.10c) 

r  =  B,  0  <  z  <  1 :  Tx  =  0 

(4.10d) 

The  transient  portion  of  the  problem  is: 

d%  1  dT2  d%  dT2 
dr2  r  dr  dz2  dt 

(4.11) 

dT 

0  <r<B,  z  =  1 :  2  =  0 

dz 

(4.11a) 

0 <r<B,  z  =  0  :  T2=  0 

(4.11b) 

dT 

r  =  0,  0  <  z  <  1 :  2  =0 

dr 

(4.11c) 

r  =  B,  0  <  z  <  1 :  T2=  0 

(4.  lid) 

0  <r<B,  0<z<l,  t  =  0:  T2  =  \-Tx(r,  z) 

(4.  lie) 
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The  steady  state  portion  of  the  problem  is  solved  by  assuming  the  functional  form 
of  the  solution  is  7j  (r,  z)  =  P(r)S(z) .  Separation  of  variables  produces  the  following 
equations  and  boundary  conditions. 


rf2?-Sr  =  0 
dz2 

(4.12) 

S(B)  =  0 

(4.12a) 

d2P  \dP  2  n 

— -  + - +  Pr  =  0 

dr  r  dr 

(4.13) 

P(B)  =  0 

(4.13a) 

®  =0 

Sr  ,.o 

(4.13b) 

The  solution  of  the  axial  direction  problem  is 

S(z)  =  a0  sinh  ruz 

(4.14) 

and  that  of  the  radial  direction  is 

P(r)  =  J0(rar) 

(4.15) 

where  the  eigenvalues  (ru's)  are  the  positive  roots  of  J0(ruB)  =  0 . 

Assembling  to  obtain 

the  solution  to  the  steady-state  problem  Tx(r,z)  =  P(r)S(z) 

and  applying  the 

nonhomogeneous  boundary  condition, 

8T. 

=~q0 

dz  Z=1 

(4.10a) 

the  following  equation  is  obtained. 

00 

-^r"  =^<,ruJ0(z‘ur)coshru 

U=1 

(4.16) 

By  using  the  principle  of  orthogonality,  the  constants  can  be  found. 

d  =  ~2q° 

Br2uJl(ruB)coshru 

(4.17) 

Thus,  the  steady-state  portion  of  the  problem  is 

T  r  _x  -‘in  ^  J0(rur)sinhTuz 

B  r\j t  (tuB)  cosh  ru 

(4.18) 

where  the  ru's  are  given  by  the  positive  roots  of  J0(tuB)  =  0 . 
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The  transient  problem  is  solved  in  a  similar  manner  by  assuming  that  the  solution 


can  be  obtained  by  separation  of  variables  by  letting  T2(r,  z,  t)  =  R(r)Z(z)T (t) .  The 
transient  dimensionless  conservation  of  energy  equation  and  the  boundary  conditions  and 
initial  condition  can  be  solved  as  three  separate  problems: 


o 

II 

l-H 

(4.19) 

d  Z  +Zrt2  =0 
dz 2 

(4.20) 

Z(  0)  =  0 

(4.20a) 

^  =0 

dz  Z=1 

(4.20b) 

d2R  1  dR  do2  a 

— r  + - +  RJ3  =  0 

dr  r  dr 

(4.21) 

R(B)  =  0 

(4.21a) 

"  =0 

dr  r=Q 

(4.21b) 

where  the  separation  constants  are  related  by  A2  =  J32  +  rf .  The 

solution  of  the 

transient  portion  of  the  problem  is 

r(0  =  Cj  exp[-(/?2  +?72)t] 

(4.22) 

The  solution  of  the  axial  direction  problem  is  given  by 

Z(z)  =  bx  sin^nz,  7n=(2n-l)y,  n  =  1,  2,  3... 

(4.23) 

The  solution  of  the  radial  direction  problem  is 

R(r)  =  J0(Pmr) 

(4.24) 

where  the  eigenvalues  \Pm's)  are  the  positive  roots  of  J0(/3mB)  =  0.  Assembling  for  the 
solution  to  the  transient  problem  T2(r,  z,  t)  =  R(r)Z(z)T (t)  and  applying  initial  condition 
T2 (r , z,0)  =  \  —  Tl(r,z)  gives 


1  ~Tx(r,z)  =  (4.25) 

m=l  n=l 

The  constants  can  be  detennined  by  using  orthogonality. 
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and  the  J3J s are  given  by  the  positive  roots  of  J0(j3mB )  =  0 .  Assembling  for  the  original 

problem  gives  the  dimensionless  temperature  distribution  within  the  right  circular 
cylinder 


PL+ni 


>/o(/?n/*)sin/7„^*-exp[-(/?2  +  r/2)t*] 


2q'o  *y  Jo(0mr*)sinhpmz* 
B  hpiJx(PmB)  cosh/?m 


(4.28) 

where  the  asterisks  have  been  re-introduced  into  the  equation  for  clarity. 

A  numerical  model  of  the  right-circular  cylinder  was  developed  using  ANSYS  5.4 
finite  element  software  for  comparison  with  the  analytical  solution.  The  thermal  model 
consisted  of  8-node,  Plane77,  axisymmetric,  quadrilateral  elements.  A  grid  independence 
study  was  performed  to  detennine  the  proper  number  of  elements  for  the  steady  state 
numerical  model.  The  uniformly-meshed  model  was  tested  for  2  x  2,  4  x  4,  and  10x10 
elements.  Figure  18  shows  the  dimensionless  temperature  along  the  radius  of  the  model 
at  three  different  z*-locations.  The  results  of  the  numerical  model  closely  match  those  of 
the  analytical  solution  for  10  x  10  elements.  The  maximum  percent  difference  between 
the  10  x  10  numerical  model  and  the  analytical  solution  was  0.05%.  A  time-step 
convergence  study  was  then  performed  on  the  transient  model  to  determined  appropriate 
values  for  At*.  Figure  19  shows  the  dimensionless  temperature  at  three  z*-locations 
along  the  centerline  of  the  right  circular  cylinder  r*  =  0.  Here  the  transient  numerical 
model  is  shown  for  time-step  values  of  At*  =  0.005,  0.01,  and  0.05.  For  At*  =  0.005,  the 
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maximum  percent  difference  between  the  numerical  and  analytical  solutions  occurred  at 
t*  =  0.005  and  (r*,  z*)  =  (0,  1),  with  a  value  of  2.69%.  The  maximum  percent  difference 
between  the  numerical  and  analytical  solutions  for  all  other  time  step  values  was  1.49%. 

Grid  Independence  Study 

A  steady  state,  numerical  model  of  the  cylindrical  glass  heater  pedestal  assembly 
was  built  in  ANSYS,  using  the  same  preferences  and  elements  described  above.  Instead 
of  a  single  thickness,  however,  the  model  consisted  of  several  layers,  reflecting  the 
structure  of  the  heater,  substrate,  and  post,  as  well  as  thermocouple  locations  within  the 
assembly.  The  height  of  the  pedestal,  excluding  the  heater  and  substrate,  was  assumed  to 
be  (//3  +  H4+  Hs)  =  10  cm  for  modeling  purposes.  The  height  of  the  glass  base  was  H5 

=  9.8  cm,  and  as  stated  in  the  Experimental  Design  and  Testing  section,  the  thickness  of 
each  borosilicate  layer  was  =  H4  =  1  mm.  The  thickness  of  the  Coming  glass 
substrate  was  H2  =  0.5  mm,  and  the  ITO  heater  was  H\  =  1  pm  thick.  A  thermal 
conductivity,  k  =  1.04  W/m-K  was  assumed  for  the  entire  model.  To  reflect  the 
symmetry  of  the  model,  a  zero-heat  flux  boundary  condition  was  placed  on  the  centerline 
of  the  structure.  The  spray  convection  was  modeled  as  a  convective  heat  transfer 
coefficient  (hi)  on  the  top  edge,  and  the  convection  from  flow  down  the  side  walls  of  the 
pedestal  assembly  was  represented  by  a  separate  convective  heat  transfer  coefficient 
along  the  right  edge  of  the  model  ( h 2).  The  bottom  of  the  pedestal  was  assumed  to  be 
adiabatic.  The  heater  was  modeled  with  a  volumetric  heat  generation  applied  to  the  top¬ 
most  layer  (ITO)  of  the  model. 

A  grid  independence  study  was  performed  to  ensure  adequate  meshing  of  the 
model.  The  coarse  mesh  was  determined  by  thermocouple  location.  In  the  coarse  mesh, 
the  maximum  element  width  in  the  r-direction  was  A r  =  0.00235  m.  In  the  z-direction,  the 
heater,  substrate,  and  the  two  thin  layers  were  each  one  element  thick,  and  the  remaining 
portion  of  the  base  was  divided  into  10  elements,  with  a  spacing  ratio  of  30,  so  that  the 
element  size  increased  downward.  This  produced  56  total  elements  in  the  coarse  model. 
The  mesh  was  refined  by  doubling  the  number  of  elements  in  both  directions,  with  the 
exception  of  the  heater,  where  its  thickness  remained  at  one  element.  The  result  was  a 
fine  mesh  consisting  of  216  elements.  Figure  20  and  Figure  21  illustrate  the  coarse  and 
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fine  meshes,  respectively.  For  the  grid  independence  study,  convective  heat  transfer 
coefficient  values  of  h\  =  10,000  W/(m-K)  and  hi  =  400  W/(m-K)  were  chosen,  along 
with  a  heater  power  of  Q  =  20  W.  Temperatures  were  recorded  at  all  nodes  and 
compared.  The  maximum  percent  difference  between  the  coarse  and  fine  meshes  was 
0.058%.  Figure  22  shows  the  temperature  along  the  radial  direction  for  three  separate  z- 
locations.  As  expected,  the  temperature  within  the  pedestal  assembly  decreases  in  the 
radial  direction  due  to  the  convection  at  the  outer  radius.  In  addition,  the  temperature 
increases  in  the  axial  direction  due  to  the  presence  of  the  ITO  heater  at  the  upper  surface. 
Based  on  these  results,  it  was  detennined  that  the  model  was  sufficiently  grid- 
independent  when  the  fine  mesh  is  used. 
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Figure  16:  Transient  heat  conduction  in  a  right  circular  cylinder:  (a)  Geometrical 
nomenclature;  ( b )  Solution  domain. 
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Figure  18:  Dimensionless  temperature  versus  dimensionless  radius  for  the  right  circular 
cylinder  (steady  state,  B  =  1,  q'0*  =  1 ). 
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Figure  19:  Dimensionless  temperature  versus  dimensionless  time  for  the  right  circular 
cylinder  (B  =  1,  r*  =  0,^*  =  1). 
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Figure  20:  Coarse  mesh  for  the  glass  heater  pedestal  assembly. 
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Figure  22:  Temperature  versus  radial  location  for  z  =  0.1,  0.099,  and  0.098  m  for  the 
coarse  and  fine  meshes. 
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V.  Results  and  Discussion 

Baseline  experiments  using  subcooled  FC-72  fluid  were  carried  out  to  determine 
the  effects  of  the  volumetric  flow  rate  of  the  spray,  the  heat  transfer  rate,  and  heater 
orientation  with  respect  to  terrestrial  gravity  on  the  heat  transfer  coefficient  on  the  top 
surface  of  a  glass  heater  pedestal  assembly.  The  ranges  of  the  parameters  in  the 
experiments  are  as  follows:  4.0  <  V  <  10.0  gal/hr,  16.0  <  Q  <  50.0  W.  Experiments  were 
executed  for  two  heater  orientations:  the  bottom  nozzle  spraying  downward  onto  an 
upward-facing  heater,  and  the  top  nozzle  spraying  upward  onto  a  downward-facing 
heater.  Raw  data  from  a  typical  experiment  is  shown  in  Figure  23.  As  can  be  seen,  the 
temperatures  within  the  pedestal  assembly  closely  follow  the  heat  input  to  the  ITO  heater 
with  respect  to  time.  The  average  volumetric  flow  rate  of  the  spray  was  relatively 
constant,  but  some  variation  around  the  average  can  be  seen,  which  was  due  to 
fluctuations  within  the  fluid  collection  system.  While  the  vapor  pressure  within  the  spray 
chamber  slowly  increased  over  the  period  of  the  test  due  to  the  heat  input  to  the  working 
fluid  through  the  glass  heater  pedestal,  the  temperature  of  the  fluid  remained  below  the 
saturation  temperature,  Tsat. 

In  order  to  determine  the  heat  transfer  coefficient  of  the  spray  on  the  ITO  heater 
surface,  the  following  procedure  was  used.  For  a  given  experimental  run  in  which  the 
volumetric  flow  rate  of  the  spray  and  the  heater  power  were  specified,  the  average 
steady-state  temperatures  for  thermocouples  1,  2  and  3  (TCI,  TC2,  TC3)  were  recorded. 
The  numerical  model  of  the  glass  heater  pedestal  assembly  was  then  used  to  match  the 
experimental  temperatures  at  the  specified  geometric  locations  of  the  thennocouples. 
This  involved  varying  both  the  spray  heat  transfer  coefficient  at  the  top  surface  of  the 
ITO  heater  (hi),  and  the  heat  transfer  coefficient  along  the  outer  radius  of  the  pedestal 
assembly  (hi).  Results  from  a  typical  comparison  between  the  experimental  data  and  the 
numerical  model  are  shown  in  Figure  24.  The  spray  heat  transfer  coefficient  was  varied 
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over  a  large  range,  while  the  heat  transfer  coefficient  on  the  outer  radius  of  the  pedestal 
assembly  was  fixed  at  a  particular  value.  The  points  at  which  the  numerical  temperatures 
at  the  locations  of  thennocouples  1,  2,  and  3  matched  those  of  the  experimental  values 
was  noted  for  that  particular  value  of  the  outer  radius  heat  transfer  coefficient.  The 
process  was  then  repeated  with  a  different  value  of  the  outer  radius  heat  transfer 
coefficient,  which  was  varied  from  0  <  h2  <  400  W/(m“-K).  This  range  was  chosen  based 
on  experimental  data  reported  in  the  open  literature  (Incropera  and  DeWitt,  2001),  and  is 
considered  to  be  a  conservative  range,  meaning  that  it  is  believed  that  the  actual  average 
heat  transfer  coefficient  is  within  this  range.  It  should  be  emphasized  that  the  numerical 
model  imposed  a  single  value  for  the  heat  transfer  coefficient  at  the  outer  radius.  In 
reality,  the  heat  transfer  is  probably  high  very  close  to  the  top  surface,  and  quickly 
decreases  to  zero  in  the  axial  direction,  and  the  heat  transfer  coefficient  varies  along  the 
wall.  Therefore,  the  assumption  that  the  range  of  the  average  heat  transfer  coefficient  lies 
between  0  <  h2  <  400  W/(m2-K)  is  justified. 

Once  the  numerical  model  was  run  for  the  entire  range  of  /i2,  the  range  of  the 
spray  heat  transfer  coefficient  was  defined  for  each  thennocouple  by  the  values  at  which 
the  numerical  temperature  values  matched  the  experimental  value.  The  range  was 
contained  between  the  intersections,  at  I12  =  0  and  hi  =  400  W/(nr-K),  of  the  numerical 
data  with  the  experimental  data.  This  is  illustrated  in  Figure  24,  with  vertical  lines  to 
show  the  range  for  each  thermocouple.  It  should  be  noted  that  the  actual  value  of  the 
spray  heat  transfer  coefficient  is  not  known,  but  it  is  believed  to  lie  in  the  range  discussed 
above.  To  simplify  the  presentation  of  the  data,  the  spray  heat  transfer  coefficient  was 
defined  as  the  average  of  the  range,  as  shown  in  Figure  25  and  Figure  26.  Here,  error  bars 
are  used  to  denote  the  extent  of  the  range.  As  can  be  seen,  the  spray  heat  transfer 
coefficient  was  not  a  strong  function  of  the  heat  input,  due  to  the  fact  that  the  fluid  was 
subcooled.  Allowing  the  fluid  temperatures  to  reach  those  greater  than  the  saturation 
temperature  would  have  created  a  greater  dependence  on  the  heat  input.  The  spray  heat 
transfer  coefficient  was  more  dependent  on  the  volumetric  flow  rate  of  the  subcooled 
fluid,  and  increased  with  the  flow  rate  of  the  spray,  as  expected.  Also  worth  noting  is  that 
Thermocouple  1  is  the  most  accurate  predictor  of  spray  heat  transfer  coefficient,  since  it 
was  located  closest  to  the  heater  surface  and  was  least  affected  by  the  wall  heat  transfer 
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coefficient,  hi,  as  demonstrated  graphically  in  Figure  24.  For  all  test  cases,  the  estimated 
spray  heat  transfer  coefficient  was  10,000  <hx<  15,000  W/(m“-K),  with  an  approximate 
rms  error  of  E(hxA )  ~  20% ,  as  discussed  in  Baysinger,  et  al.  (2004). 

It  is  worth  noting  that  the  temperatures  measured  in  the  downward-facing  heater 
pedestal  assembly  were  generally  slightly  lower  than  those  measured  in  the  upward- 
facing  heater  pedestal.  Figure  27  shows  the  behavior  of  the  spray  for  both  the  downward¬ 
facing  heater  pedestal  and  the  upward-facing  heater  pedestal.  In  both  cases  a  stagnation 
line  is  visible  within  the  sump  cap.  This  line  is  created  by  droplets  of  fluid  which  deflect 
from  the  heater  surface  to  the  inside  wall  of  the  cap.  For  the  upward-facing  heater 
pedestal,  the  stagnation  line  is  visible  slightly  above  the  heater  surface.  For  the 
downward-facing  heater  pedestal,  however,  the  stagnation  line  forms  further  up  on  the 
pedestal,  in  the  vicinity  of  the  thermocouples.  Also  visible  in  the  region  of  the  stagnation 
line  is  an  area  of  turbulent  flow.  It  is  believed  that  the  fluid  in  this  area  provides  extra 
cooling  to  the  pedestal,  resulting  in  the  slightly  lower  temperatures  observed  for  the 
downward-facing  heater  pedestal  assembly.  In  turn,  the  lower  temperatures  in  the 
regions  of  thermocouples  2  and  3  (TC2,  TC3)  result  in  the  higher  estimated  spray  heat 
transfer  coefficients  for  those  locations,  as  illustrated  in  Figure  26. 

While  most  of  the  heat  generated  was  removed  by  the  spray,  a  small  portion  was 
lost  down  the  pedestal  assembly  via  conduction.  It  was  desired  to  estimate  this  heat 
conducted  into  the  pedestal  in  order  to  validate  the  use  of  a  simplified  estimation  of  the 
heat  transfer  coefficient  as  outlined  by  Baysinger  et  al.  (2004).  The  estimated  spray  heat 
transfer  coefficient  was  used  to  determine  the  heat  lost  through  the  pedestal  due  to 
conduction,  A Q ,  for  several  experimental  tests  using  the  upward-facing  heater  pedestal. 
First,  the  actual  amount  of  heat  removed  by  the  sprayed  fluid,  Qx ,  was  found  using  the 
convection  equation: 

Qx=hl(7rRp2)(Ts(r)-Ttl), 

where  hi  is  the  estimated  spray  heat  transfer  coefficient,  Rp  is  the  radius  of  the  glass 
heater  pedestal  assembly,  Ts(r)  is  the  local  surface  temperature  given  by  the  numerical 
model,  and  Tfl  is  the  experimentally  measured  film  temperature  across  the  top  of  the 
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heater  as  described  in  the  Mathematical  Modeling  section.  The  result  was  then  used  in 
the  following  equation  to  find  A Q : 

AQ  =  \-Ql 

a 

where  Qh  is  the  power  generated  by  the  ITO  heater. 

Figure  28  displays  temperature  versus  power  for  the  average  steady  state  surface 
temperature,  Ts,  across  the  heater  and  the  temperature  of  the  surface  located  directly 
above  Thermocouple  1 ,  Ts  TC1 .  As  expected,  increasing  the  power  supplied  to  the  heater 

results  in  increased  temperature.  Increasing  the  volumetric  flow  rate  provides  more 
coolant  to  the  heated  surface,  resulting  in  decreased  temperatures  within  the  glass  heater 
pedestal  assembly,  also  shown  in  Figure  28. 

Figure  29  shows  the  percent  of  heat  lost  due  to  conduction  within  the  pedestal, 
A Q,  over  the  radius  of  the  heater  surface.  The  heat  lost  increases  approaching  the  outer 
edge,  due  to  the  convection  along  the  side  wall  of  the  pedestal  from  the  fluid  moving 
through  the  sump  structure.  Finally,  Figure  30  shows  A Q  versus  the  input  heater  power, 

based  on  Ts  and  based  on  Ts  TC1 .  As  seen  in  the  graph,  increasing  power  increases  A Q , 
while  increasing  the  volumetric  flow  rate  decreases  A Q .  In  addition,  the  heat  lost  due  to 
conduction  was  in  the  range  of  A Q  =  1  to  2.5%  which  is  comparable  to  earlier  estimates 
of  1%  (Baysinger  et  al.,  2004).  This  shows  that  the  one-dimensional  method  for 
estimating  the  spray  heat  transfer  coefficient  is  valid  as  long  as  the  thermocouple  that  is 
used  to  detennine  the  heat  transfer  coefficient  is  close  to  the  centerline  and  top  surface  of 
the  heater  pedestal.  While  this  infonnation  will  be  used  in  the  analyses  and  planning  of 
future  experimental  tests,  it  should  be  noted  that  these  estimates  are  reliable  only  for 
subcooled  conditions  of  the  fluid. 
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Figure  23:  Temperatures  within  the  upward-facing  glass  heater  pedestal  assembly,  input 
power  to  the  ITO  heater,  spray  chamber  pressure,  and  volumetric  flow  rate  of  the  spray 
versus  time  for  a  typical  test  run. 
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Figure  24:  Thermocouple  temperatures  within  the  upward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient,  h\  (V  =  5.0 gal/hr,  Q  =  20  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  25:  Spray  heat  transfer  coefficient,  hi,  versus  heater  power  (upward-facing  heater) 
(a)  V  =  5.0  gal/hr;  (b)  V  =  8.0  gal/hr;  (c)  V  =  10.0 gal/hr. 
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Figure  26:  Spray  heat  transfer  coefficient,  h\,  versus  heater  power  (downward-facing 
heater)  (a)  V  =  4.0  gal/hr;  (b)  V  =  5.0  gal/hr;  (c)  V  =  6.0  gal/hr. 
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Figure  27:  Behavior  of  the  sprayed  fluid:  (a)  Downward-facing  heater;  (b)  Upward- 
facing  heater. 
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Figure  28:  Temperature  versus  power  (a)  Average  surface  temperature;  (b)  Surface 
temperature  directly  above  thermocouple  1  (upward-facing  heater). 
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Figure  29:  Local  heat  lost  due  to  conduction  within  the  upward-facing  glass  pedestal 
heater  assembly  versus  radial  location. 
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Figure  30:  Heat  lost  due  to  conduction  in  the  upward-facing  glass  heater  pedestal 
assembly  versus  power:  (a)  Based  on  average  surface  temperature;  (b)  Based  on  surface 
temperature  directly  above  thermocouple  1 . 
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VI.  Conclusions 

Baseline  tests  were  performed  for  a  spray  cooling  system  using  subcooled  fluid 
under  terrestrial  gravity  conditions.  A  preliminary  numerical  study  was  done  to  ensure 
that  the  finite  element  software  would  provide  a  reliable  approximation  of  the 
experimental  setup.  A  steady  state  numerical  model  of  the  glass  heater  pedestal  assembly 
was  then  built  and  tested  using  the  finite  element  software.  The  spray  heat  transfer 
coefficient  along  the  top  of  the  heated  surface  was  estimated  for  several  volumetric  flow 
rates,  heat  transfer  rates,  and  heater  orientations  with  respect  to  gravity.  Several  key 
observations  were  made: 

•  The  spray  heat  transfer  coefficient,  hx ,  was  a  weak  function  of  heat  input,  Qh ,  but 

varied  significantly  with  volumetric  flow  rate,  V  . 

•  The  temperature  at  thermocouple  1  (TCI)  was  not  significantly  affected  by  h2, 
making  it  the  best  predictor  of  hx . 

•  Temperatures  within  the  heated  pedestal  remained  below  the  saturation 
temperature,  Tsat ,  in  the  spray  chamber,  indicating  that  the  assumption  of 
subcooled  fluid  was  appropriate  for  the  given  conditions. 

•  For  a  volumetric  flow  range  of  4.0<F<10.0  gal/hr  and  power  range  of 
0.0  <  Q  <  50.0  W,  hx  was  estimated  to  be  on  the  order  of  10,000  <  hx  <  15,000 
W/(m“-K)  for  both  the  upward-  and  downward-facing  heaters. 

•  Local  heat  conduction  into  the  upward-facing  heater  pedestal  assembly  was  fairly 

r 

constant  for  —  <  0.5  ,  and  was  found  to  be  1.0  <  A O  <  2.5%  . 
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Appendix  A:  Additional  Figures 


This  appendix  contains  a  complete  list  of  the  figures  used  to  determine  the  values 
presented  in  the  main  body  of  the  thesis. 
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Figure  3 1 :  Thennocouple  temperatures  within  the  upward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient:  (V  =  5.0  gal/hr,  Q  =  30  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  32:  Thermocouple  temperatures  within  the  upward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  5.0 gal/hr,  Q  =  50  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  33:  Thermocouple  temperatures  within  the  upward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  8.0 gal/hr,  Q  =  20  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  34:  Thermocouple  temperatures  within  the  upward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  8.0 gal/hr,  Q  =  30  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  35:  Thermocouple  temperatures  within  the  upward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  8.0 gal/hr,  Q  =  50  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  36:  Thermocouple  temperatures  within  the  upward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  10.0 gal/hr,  Q  =  20  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  37:  Thermocouple  temperatures  within  the  upward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  10.0 gal/hr,  Q  =  30  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 


67 


•Q—  TC 


hi  (W/m2-K) 


Figure  38:  Thermocouple  temperatures  within  the  upward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  10.0 gal/hr,  Q  =  50  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  39:  Thermocouple  temperatures  within  the  downward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  4.0 gal/hr,  Q  =  16  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  40:  Thermocouple  temperatures  within  the  downward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  4.0 gal/hr,  Q  =  30  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 


70 


hspray  (W/llf-K) 


Figure  41 :  Thermocouple  temperatures  within  the  downward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  4.0 gal/hr,  Q  =  45  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  42:  Thermocouple  temperatures  within  the  downward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  5.0 gal/hr,  Q  =  20  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  43:  Thermocouple  temperatures  within  the  downward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  5.0 gal/hr,  Q  =  30  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  44:  Thermocouple  temperatures  within  the  downward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  5.0 gal/hr,  Q  =  50  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  45:  Thermocouple  temperatures  within  the  downward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  6.0 gal/hr,  Q  =  16  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  46:  Thermocouple  temperatures  within  the  downward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  6.0 gal/hr,  Q  =  30  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Figure  47:  Thermocouple  temperatures  within  the  downward-facing  glass  heater  pedestal 
assembly  versus  spray  heat  transfer  coefficient  (V  =  6.0 gal/hr,  Q  =  45  W):  (a) 
Thermocouple  1;  (b)  Thermocouple  2;  (c)  Thermocouple  3. 
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Appendix  B:  Supplemental  Materials 


This  appendix  contains  a  listing  of  the  MATLAB  code  used  in  the  preliminary  analysis  of 
a  right  circular  cylinder,  as  well  as  ANSYS  log  and  output  fdes  of  input  and  output  used 
to  build  the  numerical  model  of  the  glass  heater  pedestal  assembly. 
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MATLAB  Program  Code 


%%Function  to  find  Tempi  for  ME890 

function  Tempi  =  Tempi 

format  long 
B=1  ; 

R=B; 

Xel=10; 

Zel=10; 

Xnodes=Xel+l ; 

Znodes=Zel+l ; 

Xinc=R/Xel ; 

Zinc=l / Zel ; 

%%Initialize  Temp2  matrix 
Temp 1  =  zeros ( Znodes+2  ,  Xnodes+2 )  ; 

%%Set  nodes  and  position  labels 
for  p=3: Xnodes+2 

Tempi ( 1 , p) =p-2 ; 

Tempi (2 , p) = (p-3 ) *Xinc; 

end 

for  q=3: Znodes+2 

Tempi (q, 1 ) =q-2 ; 

Tempi (q, 2)  =  ( q—  3 ) *Zinc; 

end 

%%Initialize  increment  matrix 
IncMat= zeros ( Znodes+2 , Xnodes+2 ) ; 

%%Set  nodes  and  position  labels 
for  p=3: Xnodes+2 

IncMat ( 1 , p) =p-2 ; 

IncMat (2 , p) = (p-3 ) *Xinc; 

end 

for  q=3: Znodes+2 

IncMat (q, 1 ) =q-2 ; 

IncMat (q, 2) = (q-3) *Zinc; 

end 

%%Get  T1  at  each  node 
for  n=3: Znodes+2 
z=Templ (n, 2 ) ; 

for  m=3: Xnodes+2 
r=Templ (2 , m) ; 


%Radial 

%Number 

%Number 

%Number 

%Number 


Dimension 
elements 
elements 
nodes  in 
nodes  in 


in  x-direction 
in  z-direction 
x-direction 
z-direction 


a— 1  ; 
b=l ; 

moveon=0 ; 
epsilon=le-6; 
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t=5  ; 

while  moveon==0 

sumold=0 ; 
oldroot=0 ; 

for  i=l:a 

root=besselroot (oldroot) ; 
beta=root/B; 

for  j=l:b 

eta= ( (2* j ) -1) * (pi/2) ; 
c— ( (1/eta)  + ( ( (-1) Aj ) 

/ (betaA2+etaA2 ) ) ) / (beta*bessel j  (1,  (beta*B) ) ) ; 

quantity= (c*besselj (0, (beta*r) ) *sin (eta*z) 
*exp (- (betaA2+etaA2 ) *t) ) ; 

sumnew=sumold+quantity; 

diff=abs ( sumnew-sumold) ; 

sumold=sumnew; 

oldroot=root; 

end 

end 

if  diff<epsilon 
moveon=l ; 

else 

a=a+l ; 
b=b+l ; 

end 

end 

IncMat (n, m) =a; 

Tempi (n, m) =sumnew*4/B; 

a; 

b; 

if  Tempi (n, m) <le-16 
Tempi  (n,  m)  =0; 

end 


end 
end 
Tempi ; 

IncMat (n,  m)  ; 
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%Function  to  find  Temp2  for  ME890 

function  Temp2  =  Temp2 

format  long 
B=1  ; 

R=B; 

Xel=10; 

Zel=10; 

Xnodes=Xel+l ; 

Znodes=Zel+l ; 

Xinc=R/Xel ; 

Zinc=l / Zel ; 

%Initialize  Temp2  matrix 
Temp2  =  zeros ( Znodes+2  ,  Xnodes+2 )  ; 

%Set  nodes  and  position  labels 
for  p=3: Xnodes+2 

Temp2 ( 1 , p) =p-2 ; 

Temp2 (2 , p) = (p-3 ) *Xinc; 

end 

for  q=3: Znodes+2 

Temp2 (q, 1 ) =q-2 ; 

Temp2 (q, 2) = (q-3) *Zinc; 

end 

%Initialize  increment  matrix 
IncMat= zeros ( Znodes+2 , Xnodes+2 ) ; 

%Set  nodes  and  position  labels 
for  p=3: Xnodes+2 

IncMat ( 1 , p) =p-2 ; 

IncMat (2 , p) = (p-3 ) *Xinc; 

end 

for  q=3: Znodes+2 

IncMat (q, 1 ) =q-2 ; 

IncMat (q, 2) = (q-3) *Zinc; 

end 

%Get  T2  at  each  node 
for  n=3: Znodes+2 
z=Temp2 (n, 2 )  ; 

for  m=3: Xnodes+2 
r=Temp2 (2 , m) ; 
epsilon=le-6; 
moveon=0 ; 
inc=l ; 
sumold=0 ; 
oldroot=0 ; 

while  moveon==0 

root=besselroot (oldroot)  ; 
beta=root/B; 


%Radial 

%Number 

%Number 

%Number 

%Number 


Dimension 
elements 
elements 
nodes  in 
nodes  in 


in  x-direction 
in  z-direction 
x-direction 
z-direction 
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sumnew=sumold+ (besselj (0, (beta*r) ) *sinh (beta*z) 

/ (betaA2*bessel j (1, (beta*B) ) *cosh (beta) ) ) 

diff=abs ( sumnew-sumold) ; 

if  diff<epsilon 
moveon=l ; 

else 

sumold=sumnew; 
oldroot=root; 
inc=inc+l ; 

end 

end 

IncMat (n, m) =inc; 

Temp2 (n, m) = (2/B) *sumnew; 

if  Temp2 (n, m) <le-10 
Temp2 (n,m) =0; 

end 

end 

end 
Temp 2 ; 
beta; 
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ANSYS  Log  File 


/BATCH 

!  /COM, ANSYS  RELEASE  5.4  UP19970828  09:47:49  10/29/2004 

!  * 

/FILNAM, StStateGlassModel 

!  * 

/NOPR 

/PMETH, OFF 
KEYW, PR_SET , 1 
KEYW, PR_STRUC, 0 
KEYW, PRJIHERM, 1 
KEYW, PR^FLUID, 0 
KEYW, PR_ELMAG, 0 
KEYW, MAGNOD, 0 
KEYW, MAGEDG, 0 
KEYW, MAGHFE, 0 
KEYW, MAGELC, 0 
KEYW, PR_MULTI, 0 
/GO 
!  * 

!  /COM, 

!  /COM, Preferences  for  GUI  filtering  have  been  set  to  display: 

!  /COM,  Thermal 

!  * 

/PREP7 

!  * 

ET, 1, PLANE 7 7 

i  * 


KEYOPT, 1,1,0 
KEYOPT, 1,3,1 
!  * 

!  * 

UIMP, 1 , EX,  ,  , 

r 

UIMP, 

1 , DENS ,  , 

r 

r 

UIMP, 

1 , ALPX,  , 

r 

r 

UIMP, 

1 , REFT,  , 

r 

r 

UIMP, 

1 ,  NUXY,  , 

r 

r 

UIMP, 

1 ,  PRXY,  , 

r 

r 

UIMP, 

1 ,  GXY ,  , 

r  r 

UIMP, 

1,MU,  ,  , 

r 

UIMP, 

1 , DAMP ,  , 

r 

r 

UIMP, 

1 ,  KXX,  , 

,i. 

04, 

UIMP, 

1  C 

r  ^  r  t  r 

r 

UIMP, 

1 ,  ENTH,  , 

r 

r 

UIMP, 

1 , HF ,  ,  , 

r 

UIMP, 

1 , EMI S ,  , 

r 

r 

UIMP, 

1 , QRATE, 

r  r 

r 

UIMP, 

1 , MURX ,  , 

r 

r 

UIMP, 

1 , MGXX ,  , 

r 

r 

UIMP, 

1 , RSVX,  , 

r 

r 

UIMP, 

1 ,  PERX,  , 

r 

r 

UIMP, 

1,VISC,  , 

r 

r 

UIMP, 

1 , SONC,  , 

r 

r 
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!  * 

K, 1,0, 0,0, 

K, 2, 0, .098,0, 

K, 3, 0, .099,0, 

K, 4 , 0 , .1,0, 

K, 5,0, .1005,0, 

K, 6, 0, . 100501, 0, 

K, 7 , .0047, .100501,0, 
K, 8 , .0047, .1005,0, 

K, 9 , .0047, .1,0, 

K, 10, .0047, .099,0, 

K, 11 , .0047, .098,0, 

K, 12, . 0047,0,0, 

K, 13, . 007,0,0, 

K, 14, .007, .098,0, 

K, 15, .007, .099,0, 

K, 16, . 007, .1,0, 

K, 17, .007, .1005, 0, 

K, 1 8 , .007, .100501,0, 
K, 19, .008, .100501,0, 
K, 20, . 008, . 1005,0, 

K, 21, . 008, .1,0, 

K, 22, . 008, .099,0, 

K, 23, . 008, .098,0, 


K, 24, . 008,0,0, 

!  klist, all, ,, coord 


LSTR, 

1, 

2 

LSTR, 

2, 

3 

LSTR, 

3, 

4 

LSTR, 

4, 

5 

LSTR, 

5, 

6 

!  LLIST, 
LSTR, 

ALL,  ,  , 

12, 

11 

LSTR, 

\ — 1 
\ — 1 

10 

LSTR, 

o 

1 — 1 

9 

LSTR, 

9, 

8 

LSTR, 

8, 

7 

!  LLIST, 
LSTR, 

ALL,  ,  , 

13, 

14 

LSTR, 

i — 1 

15 

LSTR, 

LO 
\ — 1 

16 

LSTR, 

16, 

17 

LSTR, 

17, 

18 

LSTR, 

24, 

23 

LSTR, 

23, 

22 

LSTR, 

22, 

21 

LSTR, 

21, 

20 

LSTR, 

20, 

19 

!  LLIST, 
LSTR, 

ALL,  ,  , 

6, 

7 

LSTR, 

7, 

18 

LSTR, 

CO 
\ — 1 

19 

LSTR, 

5, 

8 

LSTR, 

8, 

17 

LSTR, 

17, 

20 

LSTR, 

4, 

9 

LSTR, 

9, 

16 
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LSTR, 

\ — 1 

21 

LSTR, 

3, 

10 

LSTR, 

10, 

15 

LSTR, 

LO 

i — 1 

22 

LSTR, 

2, 

11 

LSTR, 

\ — 1 
\ — 1 

14 

LSTR, 

i — 1 

23 

LSTR, 

1, 

12 

LSTR, 

12, 

13 

LSTR, 

13, 

24 

!  LLIST, 

ALL,  ,  , 

FLST, 2, 4, 4 
FITEM, 2, 5 
FITEM, 2, 21 
FITEM, 2, 10 
FITEM, 2, 24 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 10 
FITEM, 2, 22 
FITEM, 2, 15 
FITEM, 2, 25 
AL, P51X 
!  alist,  all 
FLST, 2, 4, 4 
FITEM, 2, 15 
FITEM, 2, 23 
FITEM, 2, 20 
FITEM, 2, 26 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 4 
FITEM, 2, 24 
FITEM, 2, 9 
FITEM, 2, 27 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 9 
FITEM, 2, 25 
FITEM, 2, 14 
FITEM, 2, 28 
AL, P51X 
!  alist,  all 
FLST, 2, 4, 4 
FITEM, 2, 14 
FITEM, 2, 26 
FITEM, 2, 19 
FITEM, 2, 29 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 3 
FITEM, 2, 27 
FITEM, 2, 8 
FITEM, 2, 30 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 8 
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FITEM, 2, 28 
FITEM, 2, 13 
FITEM, 2, 31 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 13 
FITEM, 2, 29 
FITEM, 2, 18 
FITEM, 2, 32 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 2 
FITEM, 2, 30 
FITEM, 2, 7 
FITEM, 2, 33 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 7 
FITEM, 2, 31 
FITEM, 2, 12 
FITEM, 2, 34 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 12 
FITEM, 2, 32 
FITEM, 2, 17 
FITEM, 2, 35 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 1 
FITEM, 2, 33 
FITEM, 2, 6 
FITEM, 2, 36 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 6 
FITEM, 2, 11 
FITEM, 2, 34 
FITEM, 2, 37 
AL, P51X 
FLST, 2, 4, 4 
FITEM, 2, 11 
FITEM, 2, 16 
FITEM, 2, 38 
FITEM, 2, 35 
AL, P51X 
!  alist,  all 
AATT,  1,  , 

!  alist,  all 

FLST,  5,4, 4, ORDE, 4 
FITEM, 5, 5 
FITEM, 5, 10 
FITEM, 5, 15 
FITEM, 5,20 
CM, _Y, LINE 
LSEL,  ,  ,  , P51X 


CM, _Y1 , LINE 
CMSEL, ,_Y 

LESIZE , _Y1 ,  ,  ,1,1, 

CMDEL, _Y 
CMDEL, _Y1 
!  * 

FLST, 5,12, 4, ORDE, 8 
FITEM, 5,2 
FITEM, 5, -4 
FITEM, 5, 7 
FITEM, 5, -9 
FITEM, 5, 12 
FITEM, 5, -14 
FITEM, 5, 17 
FITEM, 5,-19 
CM,  Jf, LINE 
LSEL,  ,  ,  , P51X 

!  * 

CM, _Y1 , LINE 
CMSEL, ,_Y 

LESIZE, _Y1 ,  ,  ,2,1, 

CMDEL, _Y 
CMDEL, _Y1 

!  * 

FLST,  5, 4, 4, ORDE, 4 
FITEM, 5, 1 
FITEM, 5, 6 
FITEM, 5,11 
FITEM, 5,16 
CM, _Y, LINE 
LSEL,  ,  ,  , P51X 

!  * 

CM, _Y1 , LINE 
CMSEL, ,_Y 

LESIZE, _Y1 ,  ,  ,20,20, 

CMDEL, _Y 
CMDEL, _Y1 

!  * 

FLST,  5,  12, 4, ORDE, 12 
FITEM, 5,22 
FITEM, 5, -23 
FITEM, 5,25 
FITEM, 5,-26 
FITEM, 5,28 
FITEM, 5,-29 
FITEM, 5,31 
FITEM, 5, -32 
FITEM, 5,34 
FITEM, 5, -35 
FITEM, 5, 37 
FITEM, 5, -38 
CM,  Jf, LINE 
LSEL,  ,  ,  , P51X 

!  * 

CM, _Y1 , LINE 
CMSEL, ,  Y 


2,1 


LESIZE , _Y1 ,  ,  , 

CMDEL, _Y 
CMDEL, _Y1 
!  * 

FLST,  5, 6,4, ORDE, 6 
FITEM, 5,21 
FITEM, 5,24 
FITEM, 5,27 
FITEM, 5,30 
FITEM, 5,33 
FITEM, 5, 36 
CM, _Y, LINE 
LSEL,  ,  ,  , P51X 

!  * 

CM, _Y1 , LINE 
CMSEL, ,_Y 

LESIZE, _Y1 ,  ,  ,4,1, 

CMDEL, _Y 
CMDEL, _Y1 

!  * 

FLST, 5,4, 4, ORDE, 4 

FITEM, 5, 1 

FITEM, 5, 6 

FITEM, 5,11 

FITEM, 5,16 

CM, _Y1 , LINE 

LSEL,  ,  ,  , P51X 

*GET,_zl, LINE, , COUNT 

_z2  =  0 

*DO,_z5, l,_zl 
_z2=LSNEXT (_z2 ) 

*GET, _z3 , LINE, _z2 , ATTR, NDIV 
*GET, _z4 , LINE, _z2 , ATTR, SPAC 
*IF,_z3, GT, 0, THEN 
* IF, _z4 , NE , 0, THEN 
LESIZE, _z2, , ,_z3, l/_z4 
*ENDIF 
*ENDIF 
*ENDDO 
CMSEL, S , _Y1 
CMDELE, _Y1 
!  * 

/UI , MESH, OFF 
MSHAPE, 0, 2D 
MSHKEY, 1 

!  * 

FLST,  5, 3, 5, ORDE, 2 
FITEM, 5, 1 
FITEM, 5, -3 
CM, _Y, AREA 
ASEL,  ,  ,  , P51X 

CM, _Y1 , AREA 
CHKMSH, 'AREA' 

CMSEL, S,_Y 
!  * 

AMESH,  Y1 


CMDEL, _Y 
CMDEL, _Y1 
CMDEL, _Y2 

!  * 

FLST,  2, 8, 2, ORDE, 2 
FITEM, 2, 1 
FITEM, 2, -8 

BFE, P51X, HGEN,  , 1.00611ell, 1. 0061 lell, 1. 0061 lell, 1. 0061 lell, 

FLST,  5,  12, 5, ORDE, 2 
FITEM, 5, 4 
FITEM, 5, -15 
CM, _Y, AREA 
ASEL,  ,  ,  , P51X 

CM, _Y1 , AREA 
CHKMSH, 'AREA' 

CMSEL, S,_Y 
!  * 

AMESH, _Y1 

!  * 

CMDEL, _Y 
CMDEL, _Y1 
CMDEL, _Y2 

!  * 

!  elist, all, ,,0,0 

FLST, 2, 8, 4, ORDE, 4 

FITEM, 2, 1 

FITEM, 2, -5 

FITEM, 2, 36 

FITEM, 2, -38 

SFL,  P51X, HFLUX, 0,  , 

FLST,  2, 5, 4, ORDE, 2 
FITEM, 2, 16 
FITEM, 2, -20 

SFL,  P51X, CONV,  200,  ,  24.54,  , 

FLST,  2, 3, 4, ORDE, 2 
FITEM, 2, 21 
FITEM, 2, -23 

SFL, P51X, CONV, 10000,  ,23.667,  , 

!  SFLLIS,  ALL 
!  BFELIS,  ALL 

/UI , MESH, OFF 

FINISH 

/SOLU 

!  * 

ANTYPE, 0 

!  / STAT , SOLU 

SOLVE 


89 


THIS  PAGE  WAS  INTENTIONALLY  LEFT  BLANK 


90 


AFTER  YOU  HAVE  READ,  UNDERSTOOD,  AND  AGREED  TO  THE  PREVIOUS  NOTICES 
PRESS  <CR>  OR  <ENTER>  TO  CONTINUE 


*****  ANSYS  COMMAND 

LINE  ARGUMENTS 

INITIAL  JOBNAME 

=  file 

MEMORY  REQUESTED  (MB) 

=  64 

GRAPHICS  DEVICE  REQUESTED 

=  Xll 

START-UP  FILE  MODE 

=  READ 

GRAPHICAL  ENTRY 

=  YES 

LANGUAGE 

=  C 

DATABASE  SIZE  REQUESTED  (MB) 

=  16 

***  NOTE  ***  CP=  10.940  TIME= 

09:45:43 

There  are  no  parameters  and  no  abbreviations  defined. 

00022263  VERS ION=SOLARI S  RELEASE=  5.4  UP19970828 

CURRENT  JOBNAME=f ile  09:45:43  OCT  29,  2004  CP=  10.954 

/SHOW  SET  WITH  DRIVER  NAME=  Xll  ,  RASTER  MODE,  GRAPHIC  PLANES  = 

8 


RUN  SETUP  PROCEDURE  FROM  FILE=  /ansys54 /docu/start . ans 
/INPUT  FILE=  menust.tmp  LINE=  0 

/INPUT  FILE=  /ansys54 /docu/start . ans  LINE=  0 

ACTIVATING  THE  GRAPHICAL  USER  INTERFACE  (GUI).  PLEASE  WAIT... 
RESUME  ANSYS  DATA  FROM  FILE  NAME=StStateGlassModel . db 
***  ANSYS  GLOBAL  STATUS  *** 

TITLE  = 

NUMBER  OF  ELEMENT  TYPES  =  0 

0  ELEMENTS  CURRENTLY  SELECTED. 

0  NODES  CURRENTLY  SELECTED. 

ACTIVE  COORDINATE  SYSTEM 

INITIAL  JOBNAME  =  file 
CURRENT  JOBNAME  =  file 

CLEAR  ANSYS  DATABASE  AND  RESTART 

ANSYS/Mechanical/Emag  3D 

CURRENT  JOBNAME  REDEFINED  AS  StStateGlassModel 

***  NOTE  ***  CP=  33.094  TIME= 

09:48:42 

NEW  BACKUP  FILE  NAME=  StStateGlassModel . dbb . 

ALL  CURRENT  ANSYS  DATA  WRITTEN  TO  FILE  NAME=  StStateGlassModel . db 
FOR  POSSIBLE  RESUME  FROM  THIS  POINT 


MAX  ELEMENT  NUMBER  =  0 

MAX  NODE  NUMBER  =  0 

0  (CARTESIAN) 
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Preferences  for  GUI  filtering  have  been  set  to  display: 

Thermal 

*****  ANSYS  ANALYSIS  DEFINITION  ( PREP7 )  ***** 

ENTER  /SHOW, DEVICE-NAME  TO  ENABLE  GRAPHIC  DISPLAY 
ENTER  FINISH  TO  LEAVE  PREP7 

PRINTOUT  KEY  SET  TO  /GOPR  (USE  /NOPR  TO  SUPPRESS) 

ELEMENT  TYPE  1  IS  PLANE 7 7  2-D  8-NODE  THERMAL  SOLID 

KEYOPT (1-12) =  000  000  000  000 

CURRENT  NODAL  DOF  SET  IS  TEMP 
TWO-DIMENSIONAL  MODEL 

ELEMENT  TYPE  1  IS  PLANE 7 7  2-D  8-NODE  THERMAL  SOLID 

KEYOPT (1-12) =  000  000  000  000 

CURRENT  NODAL  DOF  SET  IS  TEMP 
TWO-DIMENSIONAL  MODEL 

ELEMENT  TYPE  1  IS  PLANE 7 7  AXI .  8-NODE  THERMAL  SOLID 

KEYOPT (1-12) =  001  000  000  000 

CURRENT  NODAL  DOF  SET  IS  TEMP 
AXI SYMMETRIC  MODEL 

MATERIAL  1  KXX  =  1.040000 

KEYPOINT  1  X, Y, Z  0.  0.  0.  IN 

CSYS=  0 

PLOT  KEY  POINTS  FROM  1  TO  1  BY  1 

KEYPOINT  2  X, Y, Z=  0.  0.980000E-01  0.  IN 

CSYS=  0 

PLOT  KEY  POINTS  FROM  1  TO  2  BY  1 

KEYPOINT  3  X, Y, Z=  0.  0.990000E-01  0.  IN 

CSYS=  0 

PLOT  KEY  POINTS  FROM  1  TO  3  BY  1 

KEYPOINT  4  X, Y, Z=  0.  0.100000  0.  IN 

CSYS=  0 

PLOT  KEY  POINTS  FROM  1  TO  4  BY  1 

KEYPOINT  5  X, Y, Z=  0.  0.100500  0.  IN 

CSYS=  0 

PLOT  KEY  POINTS  FROM  1  TO  5  BY  1 

KEYPOINT  6  X, Y, Z=  0.  0.100501  0.  IN 

CSYS=  0 
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PLOT  KEY 

POINTS 

FROM 

1  TO  6 

BY  1 

KEYPOINT 

CSYS=  0 

7 

X,  Y, 

Z  = 

0 . 470000E-02 

0.100501 

0. 

IN 

PLOT  KEY 

POINTS 

FROM 

1  TO  7 

BY  1 

KEYPOINT 

CSYS=  0 

8 

X,  Y, 

Z  = 

0 . 470000E-02 

0.100500 

0. 

IN 

KEYPOINT 

CSYS=  0 

9 

X,  Y, 

Z  = 

0 . 470000E-02 

0.100000 

0. 

IN 

KEYPOINT 

CSYS=  0 

10 

X,  Y, 

Z  = 

0 . 470000E-02 

0 . 990000E-01 

0. 

IN 

KEYPOINT 

CSYS=  0 

11 

X,  Y, 

Z  = 

0 . 470000E-02 

0 . 980000E-01 

0. 

IN 

KEYPOINT 

CSYS=  0 

12 

X,  Y, 

Z  = 

0 . 470000E-02 

0. 

0. 

IN 

KEYPOINT 

CSYS=  0 

13 

X,  Y, 

Z  = 

0 . 700000E-02 

0. 

0. 

IN 

PLOT  KEY 

POINTS 

FROM 

1  TO  13 

BY  1 

KEYPOINT 

CSYS=  0 

14 

X,  Y, 

Z  = 

0 . 700000E-02 

0 . 980000E-01 

0. 

IN 

KEYPOINT 

CSYS=  0 

15 

X,  Y, 

Z  = 

0 . 700000E-02 

0 . 990000E-01 

0. 

IN 

KEYPOINT 

CSYS=  0 

16 

X,  Y, 

Z  = 

0 . 700000E-02 

0.100000 

0. 

IN 

KEYPOINT 

CSYS=  0 

17 

X,  Y, 

z. 

0 . 700000E-02 

0.100500 

0. 

IN 

KEYPOINT 

CSYS=  0 

18 

X,  Y, 

z= 

0 . 700000E-02 

0.100501 

0. 

IN 

KEYPOINT 

CSYS=  0 

19 

X,  Y, 

z= 

0 . 800000E-02 

0.100501 

0. 

IN 

PLOT  KEY 

POINTS 

FROM 

1  TO  19 

BY  1 

KEYPOINT 

CSYS=  0 

20 

X,  Y, 

z. 

0 . 800000E-02 

0.100500 

0. 

IN 

KEYPOINT 

CSYS=  0 

21 

X,  Y, 

z= 

0 . 800000E-02 

0.100000 

0. 

IN 

KEYPOINT 

CSYS=  0 

22 

X,  Y, 

z= 

0 . 800000E-02 

0 . 990000E-01 

0. 

IN 
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KEYPOINT 

CSYS=  0 

23 

X,  Y, 

Z  = 

0 . 800000E-02 

0 . 980000E-01 

0. 

IN 

KEYPOINT 

CSYS=  0 

24 

X,  Y, 

z= 

0 . 800000E-02 

0. 

0. 

IN 

LINE  CONNECTS 

KEYPOINTS 

1 

2 

LINE  NO . = 

1  KP1  = 

1 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

2 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

2 

3 

LINE  NO . = 

2  KP1  = 

2 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

3 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

3 

4 

LINE  NO . = 

3  KP1  = 

3 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

4 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

4 

5 

LINE  NO . = 

4  KP1  = 

4 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

5 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

5 

6 

LINE  NO . = 

5  KP1  = 

5 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

6 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

12 

11 

LINE  NO . = 

6  KP1  = 

12 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

11 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

11 

10 

LINE  NO . = 

7  KP1  = 

11 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

10 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

10 

9 

LINE  NO . = 

8  KP1  = 

10 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

9 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

9 

8 

LINE  NO . = 

9  KP1  = 

9 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

8 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

8 

7 

LINE  NO . = 

10  KP1= 

8 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

7 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

13 

14 

LINE  NO . = 

11  KP1= 

13 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

14 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

14 

15 

LINE  NO . = 

12  KP1= 

14 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

15 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

15 

16 

LINE  NO . = 

13  KP1= 

15 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

16 

TAN  2  = 

0.0000 

1.0000 

0.0000 
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LINE  CONNECTS 

KEYPOINTS 

16 

17 

LINE  NO . = 

14  KP1= 

16 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

17 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

17 

18 

LINE  NO . = 

15  KP1= 

17 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

18 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

24 

23 

LINE  NO . = 

16  KP1= 

24 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

23 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

23 

22 

LINE  NO . = 

17  KP1= 

23 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

22 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

22 

21 

LINE  NO . = 

18  KP1= 

22 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

21 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

21 

20 

LINE  NO . = 

19  KP1= 

21 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

20 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

20 

19 

LINE  NO . = 

20  KP1= 

20 

TAN1  = 

0.0000 

-1.0000 

0.0000 

KP2  = 

19 

TAN  2  = 

0.0000 

1.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

6 

7 

LINE  NO . = 

21  KP1= 

6 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

7 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

7 

18 

LINE  NO . = 

22  KP1= 

7 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

18 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

18 

19 

LINE  NO . = 

23  KP1= 

18 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

19 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

5 

8 

LINE  NO . = 

24  KP1= 

5 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

8 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

8 

17 

LINE  NO . = 

25  KP1= 

8 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

17 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

17 

20 

LINE  NO . = 

26  KP1= 

17 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

20 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

4 

9 

LINE  NO . = 

27  KP1= 

4 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

9 

TAN  2  = 

1.0000 

0.0000 

0.0000 
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LINE  CONNECTS 

KEYPOINTS 

9 

16 

LINE  NO . = 

28  KP1= 

9 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

16 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

16 

21 

LINE  NO . = 

29  KP1= 

16 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

21 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

3 

10 

LINE  NO . = 

30  KP1= 

3 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

10 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

10 

15 

LINE  NO . = 

31  KP1= 

10 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

15 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

15 

22 

LINE  NO . = 

32  KP1= 

15 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

22 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

2 

11 

LINE  NO . = 

33  KP1= 

2 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

11 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

11 

14 

LINE  NO . = 

34  KP1= 

11 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

14 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

14 

23 

LINE  NO . = 

35  KP1= 

14 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

23 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

1 

12 

LINE  NO . = 

36  KP1= 

1 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

12 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

12 

13 

LINE  NO . = 

37  KP1= 

12 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

13 

TAN  2  = 

1.0000 

0.0000 

0.0000 

LINE  CONNECTS 

KEYPOINTS 

13 

24 

LINE  NO . = 

38  KP1= 

13 

TAN1  = 

-1.0000 

0.0000 

0.0000 

KP2  = 

24 

TAN  2  = 

1.0000 

0.0000 

0.0000 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  5) 

AREA  NUMBER  =  1 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  10) 

AREA  NUMBER  =  2 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  15) 
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AREA  NUMBER 


3 


DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  4) 

AREA  NUMBER  =  4 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  9) 

AREA  NUMBER  =  5 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  14) 

AREA  NUMBER  =  6 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  3) 

AREA  NUMBER  =  7 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  8) 

AREA  NUMBER  =  8 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  13) 

AREA  NUMBER  =  9 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  2) 

AREA  NUMBER  =  10 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  7) 

AREA  NUMBER  =  11 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  12) 

AREA  NUMBER  =  12 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  1) 

AREA  NUMBER  =  13 

DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  6) 

AREA  NUMBER  =  14 
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DEFINE  AREA  BY  LIST  OF  LINES 

(TRAVERSED  IN  SAME  DIRECTION  AS  LINE  11) 

AREA  NUMBER  =  15 

SET  ATTRIBUTES  FOR  ALL  SELECTED  AREAS 

MAT  =  1  REAL  =  0  TYPE  =  1  ESYS  =  0 

ATTRIBUTES  SET  FOR  15  AREAS  (OUT  OF  15  SELECTED) 

SET  DIVISIONS  ON  ALL  PICKED  UNMESHED  LINES 

TO  NDIV  =  1,  SPACING  RATIO  =  1.000 

SET  DIVISIONS  ON  ALL  PICKED  UNMESHED  LINES 

TO  NDIV  =  2,  SPACING  RATIO  =  1.000 

SET  DIVISIONS  ON  ALL  PICKED  UNMESHED  LINES 

TO  NDIV  =  10,  SPACING  RATIO  =  30.000 

SET  DIVISIONS  ON  ALL  PICKED  UNMESHED  LINES 

TO  NDIV  =  2,  SPACING  RATIO  =  1.000 

SET  DIVISIONS  ON  ALL  PICKED  UNMESHED  LINES 

TO  NDIV  =  4,  SPACING  RATIO  =  1.000 

PRODUCE  ALL  QUADRILATERAL  ELEMENTS  IN  2D. 

USE  THE  MAPPED  MESHER. 

GENERATE  NODES  AND  ELEMENTS  IN  ALL  PICKED  AREAS 

**  Meshing  of  area  1  in  progress  ** 

**  AREA  1  MESHED  WITH  4  QUADRILATERALS,  0  TRIANGLES  ** 

***  WARNING  ***  CP=  93.810  TIME= 

10:03:17 

Quadrilateral  element  1  has  an  aspect  ratio  of  1175,  which  exceeds  the 
warning  limit  of  20. 

***  WARNING  ***  CP=  93.867  TIME= 

10:03:17 

Quadrilateral  element  2  has  an  aspect  ratio  of  1175,  which  exceeds  the 
warning  limit  of  20. 

***  WARNING  ***  CP=  93.888  TIME= 

10:03:17 

Quadrilateral  element  3  has  an  aspect  ratio  of  1175,  which  exceeds  the 
warning  limit  of  20. 

***  WARNING  ***  CP=  93.923  TIME= 

10:03:17 

Quadrilateral  element  4  has  an  aspect  ratio  of  1175,  which  exceeds  the 
warning  limit  of  20. 

**  Meshing  of  area  1  completed  **  4  elements. 

**  Meshing  of  area  2  in  progress  ** 

**  AREA  2  MESHED  WITH  2  QUADRILATERALS,  0  TRIANGLES  ** 

***  WARNING  ***  CP=  94.164  TIME= 

10:03:17 
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Quadrilateral  element  5  has  an  aspect  ratio  of  1150,  which  exceeds  the 
warning  limit  of  20. 

***  WARNING  ***  CP=  94.214  TIME= 

10:03:17 

Quadrilateral  element  6  has  an  aspect  ratio  of  1150,  which  exceeds  the 
warning  limit  of  20. 

**  Meshing  of  area  2  completed  **  2  elements. 

**  Meshing  of  area  3  in  progress  ** 

**  AREA  3  MESHED  WITH  2  QUADRILATERALS,  0  TRIANGLES  ** 

***  WARNING  ***  CP=  94.356  TIME= 

10:03:18 

Quadrilateral  element  7  has  an  aspect  ratio  of  500,  which  exceeds  the 
warning  limit  of  20. 

***  WARNING  ***  CP=  94.400  TIME= 

10:03:18 

Quadrilateral  element  8  has  an  aspect  ratio  of  500,  which  exceeds  the 
warning  limit  of  20. 

**  Meshing  of  area  3  completed  **  2  elements. 


NUMBER  OF  AREAS  MESHED  =  3 
MAXIMUM  NODE  NUMBER  =  43 
MAXIMUM  ELEMENT  NUMBER  =  8 

PRODUCE  ELEMENT  PLOT  IN  DSYS  =  0 


SPECIFIED  BODY  FORCE  HGEN  FOR  ALL  PICKED  ELEMENTS  AT  STARTING  LOCATION 

1 

1 . 006110000E+11  1 . 006110000E+11  1 . 006110000E+11  1 . 006110000E+11 

***  NOTE  ***  CP=  100.908  TIME= 

10:05:11 

DELETED  BACKUP  FILE  NAME=  StStateGlassModel . dbb . 

***  NOTE  ***  CP=  100.915  TIME= 

10:05:11 

NEW  BACKUP  FILE  NAME=  StStateGlassModel . dbb . 

ALL  CURRENT  ANSYS  DATA  WRITTEN  TO  FILE  NAME=  StStateGlassModel . db 
FOR  POSSIBLE  RESUME  FROM  THIS  POINT 

GENERATE  NODES  AND  ELEMENTS  IN  ALL  PICKED  AREAS 

**  Meshing  of  area  4  in  progress  ** 

**  AREA  4  MESHED  WITH  8  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  4  completed  **  8  elements. 

**  Meshing  of  area  5  in  progress  ** 

**  AREA  5  MESHED  WITH  4  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  5  completed  **  4  elements. 

**  Meshing  of  area  6  in  progress  ** 

**  AREA  6  MESHED  WITH  4  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  6  completed  **  4  elements. 
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*  * 


*  * 


Meshing  of  area  7  in  progress 

**  AREA  7  MESHED  WITH  8  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  7  completed  **  8  elements. 

**  Meshing  of  area  8  in  progress  ** 

**  AREA  8  MESHED  WITH  4  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  8  completed  **  4  elements. 

**  Meshing  of  area  9  in  progress  ** 

**  AREA  9  MESHED  WITH  4  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  9  completed  **  4  elements. 

**  Meshing  of  area  10  in  progress  ** 

**  AREA  10  MESHED  WITH  8  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  10  completed  **  8  elements. 

**  Meshing  of  area  11  in  progress  ** 

**  AREA  11  MESHED  WITH  4  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  11  completed  **  4  elements. 

**  Meshing  of  area  12  in  progress  ** 

**  AREA  12  MESHED  WITH  4  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  12  completed  **  4  elements. 

**  Meshing  of  area  13  in  progress  ** 

**  AREA  13  MESHED  WITH  80  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  13  completed  **  80  elements. 

**  Meshing  of  area  14  in  progress  ** 

**  AREA  14  MESHED  WITH  40  QUADRILATERALS,  0  TRIANGLES  ** 

**  Meshing  of  area  14  completed  **  40  elements. 

**  Meshing  of  area  15  in  progress  ** 

**  AREA  15  MESHED  WITH  40  QUADRILATERALS,  0  TRIANGLES  ** 

***  WARNING  ***  CP=  105.484  TIME= 

10:05:32 

Quadrilateral  element  177  has  an  aspect  ratio  of  29.87,  which  exceeds 
the  warning  limit  of  20. 

***  WARNING  ***  CP=  105.531  TIME= 

10:05:32 

Quadrilateral  element  178  has  an  aspect  ratio  of  29.87,  which  exceeds 
the  warning  limit  of  20. 

***  WARNING  ***  CP=  105.555  TIME= 

10:05:32 

Quadrilateral  element  179  has  an  aspect  ratio  of  25.51,  which  exceeds 
the  warning  limit  of  20. 

***  WARNING  ***  CP=  105.590  TIME= 

10:05:33 

Quadrilateral  element  180  has  an  aspect  ratio  of  25.51,  which  exceeds 
the  warning  limit  of  20. 
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TIME= 


***  WARNING  ***  CP=  105.630 

10:05:33 

Quadrilateral  element  181  has  an  aspect  ratio  of  21.79,  which  exceeds 
the  warning  limit  of  20. 

* 

Error  messages  discontinued  after  5  messages  were  displayed. 

More  may  exist.  See  (  file. err  )  for  suppressed  messages. 

***icic**i?***i?ic**i?*i?-k'k-k'k'k'k-k'k'k'k-k'k'k'k'k'k'k'k'k'k'k'k-k'k'k'k'k'k'k'k'k'k-k'k'k'k'k'k'k'k'k'k'k'k'k'k-k'k'k'k-k'k'k 

* 

**  Meshing  of  area  15  completed  **  40  elements. 

NUMBER  OF  AREAS  MESHED  =  12 

MAXIMUM  NODE  NUMBER  =  719 

MAXIMUM  ELEMENT  NUMBER  =  216 

PRODUCE  ELEMENT  PLOT  IN  DSYS  =  0 

SURFACE  LOAD  ON  ALL  PICKED  LINES 
LOAD  LABEL  =  HFLU 

VALI ,  VALJ  =  0.  0. 

SURFACE  LOAD  ON  ALL  PICKED  LINES 
LOAD  LABEL  =  CONV 

VALI,  VALJ  =  200.00  200.00 

VAL2I,  VAL2J  =  24.540  24.540 

CONV  LOAD  SURFACE  DISPLAY  KEY  =  1 

SURFACE  LOAD  ON  ALL  PICKED  LINES 
LOAD  LABEL  =  CONV 

VALI,  VALJ  =  10000.  10000. 

VAL2I,  VAL2J  =  23.667  23.667 

***  NOTE  ***  CP=  122.061  TIME= 

10:09:06 

DELETED  BACKUP  FILE  NAME=  StStateGlassModel . dbb . 

***  NOTE  ***  CP=  122.068  TIME= 

10:09:06 

NEW  BACKUP  FILE  NAME=  StStateGlassModel . dbb . 

ALL  CURRENT  ANSYS  DATA  WRITTEN  TO  FILE  NAME=  StStateGlassModel . db 
FOR  POSSIBLE  RESUME  FROM  THIS  POINT 

*****  ROUTINE  COMPLETED  *****  CP  =  123.575 

*****  ANSYS  SOLUTION  ROUTINE  ***** 

PERFORM  A  STATIC  ANALYSIS 
THIS  WILL  BE  A  NEW  ANALYSIS 

*****  ANSYS  SOLVE  COMMAND  ***** 

TRANSFER  SOLID  MODEL  BOUNDARY  CONDITIONS  TO  FINITE  ELEMENT  MODEL 
SURFACE  LOADS  TRANSFERRED  FROM  LINES  =  70 
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TIME= 


***  NOTE  ***  CP=  126.899 

10:09:31 

There  is  no  title  defined  for  this  analysis. 

***  WARNING  ***  CP=  127.052  TIME= 

10:09:32 

Previous  testing  revealed  that  14  of  the  216  selected  elements  violate 
shape  warning  limits.  To  review  warning  messages,  please  see  the 
output  or  error  file,  or  issue  the  CHECK  command. 

***  NOTE  ***  CP=  127.090  TIME= 

10:09:32 

The  model  data  was  checked  and  warning  messages  were  found. 

Please  review  output  or  errors  file  (  file. err  )  for  these  warning 
messages . 

SOLUTION  OPTIONS 


PROBLEM  DIMENSIONALITY . AXISYMMETRIC 

DEGREES  OF  FREEDOM . TEMP 

ANALYSIS  TYPE . STATIC  ( STEADY-STATE ) 

***  NOTE  ***  CP=  127.325  TIME= 

10:09:33 

Axisymmetric  elements  present.  Nodal  loads  input  with  the  F  commands 
and  the  reaction  forces  are  on  a  full  circle  basis.  This  differs 
from 

pre-ANSYS  5.0  interpretation. 

***  NOTE  ***  CP=  127.332  TIME= 

10:09:33 

Present  time  0  is  less  than  or  equal  to  the  previous  time. 

Time  will  default  to  1. 


***  WARNING  ***  CP=  127.353  TIME= 

10:09:34 

Element  1  face  4  convection  is  on  or  near  a  zero  radius  face. 

***  WARNING  ***  CP=  127.401  TIME= 

10:09:34 

Element  9  face  4  convection  is  on  or  near  a  zero  radius  face. 

***  WARNING  ***  CP=  127.474  TIME= 

10:09:39 

Element  13  face  4  convection  is  on  or  near  a  zero  radius  face. 

***  WARNING  ***  CP=  127.500  TIME= 

10:09:39 

Element  25  face  4  convection  is  on  or  near  a  zero  radius  face. 


ici?icic*i?&ic*i?*ic*i?*ic*-k-k'k'k-k'k'k'k-k-k'k'k-k'k'k'k-k-k'k'k-k'k'k'k-k-k'k'k-k'k'k'k'k'k'k'k-k'k'k'k-k'k'k'k-k'k'k'k-k-k'k'k-k'k 

* 

Error  messages  discontinued  after  5  messages  were  displayed. 

More  may  exist.  See  (  file. err  )  for  suppressed  messages. 

* 
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TIME= 


***  NOTE  ***  CP=  127.748 

10:09:39 

The  step  data  was  checked  and  warning  messages  were  found. 

Please  review  output  or  errors  file  (  file. err  )  for  these  warning 
messages . 

***  NOTE  ***  CP=  127.754  TIME= 

10:09:39 

Results  printout  suppressed  for  interactive  execute. 

LOAD  STEP  OPTIONS 


LOAD  STEP  NUMBER . 

TIME  AT  END  OF  THE  LOAD  STEP . 

NUMBER  OF  SUBSTEPS . 

STEP  CHANGE  BOUNDARY  CONDITIONS  .... 
AXI SYMMETRIC  HARMONIC  LOADING  PARAMETERS 

MODE . 

I SYM . 

PRINT  OUTPUT  CONTROLS  . 

DATABASE  OUTPUT  CONTROLS . 


1 

1.0000 

1 

NO 

0 

1 

NO  PRINTOUT 

ALL  DATA  WRITTEN 

FOR  THE  LAST  SUBSTEP 


Range  of  element  maximum  matrix  coefficients  in  global  coordinates 
Maximum=  89.6753157  at  element  6. 

Minimum=  1 . 100479334E-02  at  element  113. 

***  ELEMENT  MATRIX  FORMULATION  TIMES 
TYPE  NUMBER  ENAME  TOTAL  CP  AVE  CP 

1  216  PLANE 7 7  0.938  0.004 

Time  at  end  of  element  matrix  formulation  CP=  129.610037. 

Estimated  number  of  active  DOF=  719. 

Maximum  wavefront=  54. 


Time  at  end  of  matrix  triangularization  CP=  130.78724. 
Equation  solver  maximum  pivot=  46.3889868  at  node  31  TEMP. 
Equation  solver  minimum  pivot=  7 . 0177546E-04  at  node  3  TEMP. 


***  ELEMENT  RESULT  CALCULATION  TIMES 
TYPE  NUMBER  ENAME  TOTAL  CP  AVE  CP 

1  216  PLANE 7 7  0.784  0.004 

***  NODAL  LOAD  CALCULATION  TIMES 
TYPE  NUMBER  ENAME  TOTAL  CP  AVE  CP 


CUM  ITER  =  1 

NEW  TRIANG 


***  PROBLEM  STATISTICS 

ACTUAL  NO.  OF  ACTIVE  DEGREES  OF  FREEDOM  =  719 

R.M.S.  WAVEFRONT  SIZE  =  36.8 


1  216  PLANE 7 7  0.183  0.001 

***  LOAD  STEP  1  SUBSTEP  1  COMPLETED. 

***  TIME  =  1.00000  TIME  INC  =  1.00000 

MATRIX 
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k  k  k 


ANSYS  BINARY  FILE  STATISTICS 
BUFFER  SIZE  USED=  16384 

0.188  MB  WRITTEN  ON  ELEMENT  MATRIX  FILE:  StStateGlassModel . emat 
0.062  MB  WRITTEN  ON  ELEMENT  SAVED  DATA  FILE: 

StStateGlassModel . esav 

0.250  MB  WRITTEN  ON  TRI ANGULARI ZED  MATRIX  FILE: 
StStateGlassModel . tri 

0.250  MB  WRITTEN  ON  RESULTS  FILE:  StStateGlassModel . rth 

***  NOTE  ***  CP=  133.681  TIME= 

10:10:03 

DELETED  BACKUP  FILE  NAME=  StStateGlassModel . dbb . 

***  NOTE  ***  CP=  133.687  TIME= 

10:10:04 

NEW  BACKUP  FILE  NAME=  StStateGlassModel . dbb . 

ALL  CURRENT  ANSYS  DATA  WRITTEN  TO  FILE  NAME=  StStateGlassModel . db 
FOR  POSSIBLE  RESUME  FROM  THIS  POINT 

COPY  INPUT  COMMANDS  FROM  DATABASE  TO  FILE  StStKlog 
FILTER  OPTION=  COMMENT 
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